UNIVERSITY OF HAW 
RARY ages 


PHILOSOPHICAL 
MAGAZINE 


FIRST PUBLISHED IN 1798 


L. 46 SEVENTH SERIES No. 383 December 1955 


A Journal of 


Theoretical Experimental 


and Applied Physics 


EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Sc. F.R.S. 
SIR GEORGE THOMSON, M.A., D.8c., F.R.S. 
PROFESSOR A. M. TYNDALL, ©.B.E., D.Sc., F.R.S. 


PRICE 15s. 0d. 


Annual Subscription £8 Os. 0d. payabie in advance 


) AND PUBLISHED BY TAYLOR & FRANCIS LTD., RED LION..COURT, FLEET ST., LONDON, E.C.4 


Taylor & Francis, Ltd., Red Lion Court, London, E.C.4 announce the 
publication of a new scientific journal, to appear in six parts per yearly 
volume 


Journal of Fluid Mechanics 


Editor : 


Dr. G. K. BATCHELOR, Cavendish Laboratory, University of Cambridge, 
Cambridge, England 


Associate Editors : 


Prof. G. F. CARrRtgR, Pierce Hall, Harvard University, Cambridge 38, 
Massachusetts, U.S.A. 


Prof. W. C. GRIFFITH, Palmer Physical Laboratory, Princeton 
University, Princeton, New Jersey, U.S.A. 


Prof. M. J. LIGHTHILL, Department of Mathematics, The University, 
Manchester, England 


The Journal will publish papers describing theoretical or experimental 
investigations of any aspect of fluid mechanics, either fundamental or 
arising in the context of natural and applied sciences such as aeronautics, 
astrophysics, chemical engineering, hydraulics, and meteorology. 
The Journal is intended to be devoted primarily to fluid mechanics 
per se, but papers describing engineering applications or relevant 
mathematical techniques will be accepted provided they are not too 
specialized in their appeal to readers. The editors propose to do all 
they can to encourage clarity of exposition in papers published in the 
Journal. 

Papers from all countries will be welcome, and proofs of papers 
from abroad will be despatched by airmail. Authors will be entitled 
to receive 50 off-prints of a paper in the Journal free of charge. 


6 parts per volume—£1 0s. 0d. ($3.00) per part 
Subscription price per volume £5 10s. Od. ($16.50) post free, payable in advance 


The first number of the Journal will appear early in 1956 


i 
— 
bo 
=p) 
Ou 

S| 


CXXXVITT. Non-Linear Regenerative Extraction of 
Synchrocyclotron Beams 


By K. J. Le Courrur 
Department of Theoretical Physics, Liverpool University 


and 8. Lipron 
Department of Statistics, Rothamsted Experimental Station* 


[Received September 12, 1955] 


ABSTRACT 


To extract a beam of a maximum energy from a synchrocyclotron, 
the deflector must be placed at the edge of the magnet, where the magnetic 
field falls off very rapidly and non-linearly with radius. This presents 
an analytical problem of coupled non-linear differential equations with 
periodic coefficients which is solved approximately. The results have 
been checked by electronic computation. It is concluded that a non- 
linear deflector can extract the beam with higher energy and perhaps 
greater theoretical efficiency than the linear one at present working at 
Liverpool. The new method is mechanically simpler than the old for 
it has no peeler, only a regenerator. The mathematics has interest 
apart from this application, for the difficulty of handling these coupled 
cyclotron equations is a major obstacle to the development of new 
particle accelerators of the fixed field type. 


§ 1. INTRODUCTION 


Tue successful operation of the regenerative beam extractor in the 
Liverpool synchrocyclotron (Le Couteur 1951, 1953, Crewe and Le Couteur 
1955) showed that the behaviour of the system was predictable and stable. 
However, the design was rather cautious and the high efficiency of 
operation entailed some sacrifice of energy; particles were extracted 
from 69 in. radius although the natural blow up of the internal beam 
occurs only at 71-5in. Further analysis suggests that this limitation 
can be overcome but the system is so strongly non-linear that precise 
predictions require extensive numerical work. 


§ 2. Equations oF MoTIon 


Let @ denote the azimuth and p and z the radial and vertical 
displacements of a particle from a synchronous orbit of radius 79, where 
the magnetic field is Hy. At the position of the particle the vertical 


* Communicated by the Authors. 
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magnetic field is H=H(r, 2) and, if the conservation of energy is used 


to express d?/dé? as 
do\? ldrd 
ay pate Site ee 
d?e= (3) (a /d6 whey ii) 


the radial Lorentz equation of motion becomes 


dp | r(rH—roH,) 1 {fdp\?__ (dz\? 
dit rH, on do) a6) 


with neglect of a term of third order in the velocities. This may be 
vritten 

dp ro 2p p°H 1 dp\? dz\2 
—— _—— = — = — — — 2 
age Pt H)( 14+ 4 + Hy 2ry \\ab) — \a0 2) 


or, if velocity 2 terms are neglected and H expanded by Taylor’s theorem, 
using d?H/dr?4+-d*H/dz*=0, 


dp r 0H 2p aeons a | rer fir! 
Ta tet (a =) e(1+ =) (a Ga) Oe hig aimee (3) 


The vertical equation of motion is 
d*z Ildpdz_ rH" 


ii db dha Hs ie te (4) 
or, neglecting velocity ? terms and using 0H"/dz=0H7?/dr, 
2 2 
Roe e+) = 0.11 ty eee Gene 
Further expansion yields the second order. equations 
dz teOn, reH « 20H 
m7 5) PH Bry: sh =o. ° . . . (6) 


These forms of the equations of motion show that the radial motion 
depends on (H—H,)/p and the vertical on dH/dr ; these quantities are 
significantly different in the present application. 

Figure 1 shows a graph of the field H at the synchrocyclotron median 
plane and of n=—(r/H)(0H/er) and Hy. Let us consider deflector 
fields which start to act at 7y>=71in., so allowing normal acceleration 
to continue almost to the natural blow-up point. The analytical problem 
is clarified if the foregoing equations of motion are written explicitly as 

d? a2 
sa + 2%=0, ae +w%=0 for p<0, . . . (7) 


; dz 
age + 2*p(1—0-152p)+0-132=0, Fa + (w?+-0-26p)z=0 


for p> O55 we a) 
where {2?=0:955=1—n and w?=0-045=n define the frequencies of 
radial and vertical oscillation in the linear region p<0. 

On the median plane the radial restoring force vanishes when 0-152p=1 
or p=6-6in. This can be seen directly from fig. 1, for the values of 
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Hy at 71 in. and 77-6 in. are equal, so a particle with energy appropriate 
to 71 in. can also circulate on a metastable orbit at 77-6 in. and would 
have to be deflected beyond this radius in order to spiral out from the 
magnet without a magnetic channel. Such a large deflection does not 
seem to be advisable, so the problem is to swing the particles into a 
magnetic channel placed at about 74 in. 


Fig. 1 
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69 7\ as) Wise 2G We 
Radial variation of H, Hr and n=—(r/H)(0H/¢r). 


The deflector field H ,, of angular width @, at the regenerator vanishes 
at 7) and increases outwards. According to (2) and (5) the regenerator 
produces impulses 


dp +) ro ( p= 
— —(—) =—6,— 47(14+—) =—-T,p, 
an i in tal 0 : "9 Z 


dz dz No ae E) a 
= = a a ere, 1 — 2= fe t a , 9) 
ai. (3). Y - Jet. or ae To ( 


and a peeler, if used, acts similarly with S replacing ?. Evidently 7’. 
and 7’, are unequal if H varies nonlinearly with p ; it will be shown that 
this feature makes it possible to combine vertical stability with large 
radial deflections. 
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The equations formulated above are too complicated for exact solution 
and an adequate numerical survey of the problem would be a very large 
task, even for an electronic computer. So an approximate analysis of the 
strongly non-linear case has been used to determine suitable deflector 
fields and the expected performance was checked by a reasonable amount 
of orbit computation. 


§ 3. RapiaL Motion 

The plot of H against r is approximated by two straight lines joining 
at yp—=41 in; 

Then the radial motion of a particle with p<0 is a simple harmonic 
motion with frequency Q, defined by (7), and for p>0, is a simple harmonic 
motion with a lower frequency 2’. Encounters with the regenerator 
all take place with p>0. 

Consider the radial motion of a particle which crosses the radius 7 
at azimuth g after the regenerator, moving inwards as in fig. 2. After 


Fig. 2 
REGENERATOR 


k=77r/N 


The dashed line indicates the path of a particle deflected by the regenerator 
from the synchronous orbit of radius 7. 


an angle k=r/Q-=7 the particle again crosses rg moving outwards with 
its original radial velocity reversed. This is exactly the effect which 
would be produced by traversing an angle k’=7/Q’ through a region 
of characteristic radial frequency 2’. After a further angle h the particle 
re-enters the regenerator. Thus the regenerative radial motion for 
p>9 may be calculated using the radial frequency ’ everywhere, if the 
gap f=g+n/Q2-+-h is replaced by 
f'=gt7/Q'+h=ft+n]/Q’—n/2. . . . . . (10) 
Accordingly the radial gain factor e4 may be calculated from (14) 
of 1951. Results are shown in fig. 3 for n=0-0625, Q=0-9683 and various 
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values of 2’. The broken curve shows that a peeler of strength S=0-2 
placed 60° before the regenerator, does not much increase the radial 
gain and further calculation shows that it does not improve the vertical 


Ome. 38 0490 590.68707 7087709 1-077 


Radial gain ¢4 as a function of regenerator strength 7’, for Q—0-9683 and 
several values of Q’. The broken curve shows effect of a peeler placed 
60° before the regenerator. 


focussing either. So it is best to work with a regenerator only, relying 
on the natural decrease of H in the region outside ry for the peeler effect. 
In this case formulae become very simple, (14) of 1951 reduces to 


fi 
Len! Opes ec LPs ee ae 
cosh A=cos {2'f sp sin COM a (11) 
Fig. 4 
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Phase angle g with conditions as in fig. 3. 


with f'=20-+n]/ Q’—n/Q. Tf the azimuth @ is measured from the centre 
of the regenerator, after capture in the deflector the motion of a particle 
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leaving the regenerator for the mth time is 


p=ce™4 sin (¢+2'0) for p>0, . . . - - (12) 
with ¢ determined from 1953, eqn. (7), and 1951, eqn. (15) as 
cot d=(e4— cos Q'f)/sin Qf’. .-. . . ~ (18) 


So in this motion p vanishes at azimuth g=(z—¢)/' and later at g--7/2. 
Afterwards with p positive the particle approaches the regenerator again 
on an orbit 

p=c’e™ gin Q'(6—g—m/@2). . . .. ss tae) 
Values of g are plotted in fig. 4. 


§ 4. VerticaL Motion 

Careful consideration of the vertical motion is necessary to avoid 
vertical instability of the beam and for high efficiency the vertical 
expansion must be minimized. 

Initially, when the amplitude of radial oscillation remains small, 
the third term in (6) may be neglected and the 1951 theory of constant 
n applies. Later when the particle is captured into the regenerative 
motion (12), (14), p is known as a function of @ so (5) or (6) become 
essentially linear differential equations with periodic coefficients. Then, 
in the approximation of § 3, the frequencies of vertical oscillation are 


w=1/n for p<0O im region k, 
(15) 


, 


w for p>0O in regions g and h. 


In the absence of deflector fields the vertical motion of a particle 
constrained to follow the radial path (5), (6), becomes unstable, because 
there is excessive disparity between w’ and w when the amplitude of 
radial motion is large enough to make w’ >0-55. However, the regenerator 
can be designed to maintain vertical stability by compensating for the 
large decrease of H in the regions g and h. 

As in (1951) eqn. (7) or Appendix A the effect on the vertical motion 
of one revolution through the deflector is given by 

“ F | (16) 


z 2Z . 
(«itt nas 4 (cal with A—GQHK= |" 


where z,, is the vertical displacement at the mth passage through the 
radial node at the meet of gy and k and G, H, K are matrices like D, F 
of 1951. 

Vertical stability requires A to have eigenvalues of the form e+*4, so 
cos A=3(~-+5) 


: hy ’ ’ (w’—w)? : : 
=c0s (wk+w’g+oa'h)— Joa? SB wk sin w'(g+h) 
ae fi . k , '} (w’—w)? : ’ 
3,7) (wk-+w'g+o'h)+ ale sin wk cos w'(g+-h) 


w2—w? 2 A 
+ Sear sin wk cos w ia Sess ly) 


Regenerative Extraction of Synchrocyclotron Beams 1271 


should be numerically less than unity. When this condition is satisfied 
the deflector increases the initial vertical height h, of the beam by a factor 


h 


max 


i =(1-++cos e/1—cos e)"?. naire > (18) 
0 


The vertical defocussing parameter cos « must be kept small and, for the 
motion in gap k is calculated, as in Appendix A of 1951, from 


and similar formulae apply to gaps g and / with w’ replacing w; cos «=1 

corresponds to the limits cos A=-+41 of vertical stability. The values of 

f and y defined by (16) give 

(w'—w)? 
2aw’ 

(w’—w)? 


$(Bo—y/w)=sin (wk+o'g+o'h)+ cos wk sin w'(g+h) 


cos (wk+w'g+u’h)+ cos wk cos w'(g+h) 


T 2wo' | 
Qu’ Tg pean: 
| ear 


+ Joa’? 8 wk cos w'(h—-g). 


19) 


O 0.5 1.O Pon af 


Vertical defocussing, cos «, for w=0-25 and various values of w’ and g. 


Putting 7’=0 in eqn. (17) shows that, in the absence of a regenerator, 
for w==0-25 vertical instability sets in when w’ > 0-55 so making cos A< —I. 
However the values of cos « plotted in fig. 5 show that vertical stability, 
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and even small vertical defocussing, can be maintained by proper choice 
of T. The values of cos e’ are similar, but rather less as the vertical 
focussing forces are stronger in gaps g and / than in k. 

Similar analysis shows that in the difficult region of large w’ vertical 
defocussing is not reduced by use of a second deflector field, whether 
placed before or after the regenerator. So these complications are not 
discussed in detail. 


§ 5. FREQUENCIES OF RADIAL AND VERTICAL MOTION 


To make the best use of the results of § 4 the effective frequencies £2’ 
and w’ of radial and vertical motion in the region p>0 must be carefully 
chosen. 

For p>0 and small z, the quarter period of the radial motion has been 
- determined from (8) by means of an elliptic integral and, expressed as 
77/222’, determines the effective frequency {’ as a function of the amplitude 
a of outward radial oscillation. These results are graphed in fig. 6. 


Fig. 6 
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Effective frequencies Q’ and a’ of radial and vertical oscillation during an 


outward swing of amplitude a. 


It does not seem necessary to give details of integration since the intuitive 
approximation 


quarter yeriod = — : aeteekeo te ey 2 
a # QT Qo 1622) Depew ihe): 


is almost as good. 
Equation (8) determines the average vertical oscillation frequency 
during an outward radial half cycle of amplitude a as approximately 
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1 TE 
-| (w?+-0-26a sin 6)?d9=0-24w+0-76(w2+0-26a)¥?2, . (21) 
0 


which is also plotted in fig. 6. 


§ 6. CHOICE OF DEFLECTOR FIELDS 
It is now possible to choose a suitable deflector field by combining the 
previous results. 
At zero radial amplitude fig. 3 shows that 7’,.>0-2 is required to initiate 
the regenerative motion. 
At 2-5in. radial amplitude, where w’=0-7, fig. 5 shows that to limit 
the vertical defocussing regenerator strength 7',>1 and phase g>60° 


are desirable, at this radius 2’=—0-8 and fig. 4 shows that the condition 
on g requires 7',<0-7. 


Fig. 7 


O naw ie 
Amplitude a 


Radial gain e4, phase g and vertical defocussing cos « as functions of the amplitude 
a of outward radial oscillation for two regenerator fields. (a), full line 
T—V—0-2. (bd), broken line 7’=0-2. V=0:25. 


These various requirements are compatible. onsider a regenerator 
with 


opp H a= Tet Ver: eet (22) 
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according to definitions (9), then 
2T 2T SS 
0 0 


The behaviour of a system with 7’=V=0-2 is shown in fig. 7(@), which 
indicates how radial gain e4, phase g and vertical defocussing cos « vary 
with the amplitude of the outward radial motion. When the amplitude 
reaches 2-5in. the radial gain is sufficient for entry into a thin-walled 
magnetic channel and the vertical defocussing is small enough for a high — 
efficiency. After an initial decrease, the phase g does not vary much 
which implies that particles entering the channel should appear to 
diverge from a point source at azimuth g+-7/ 2 and incidentally validates 
the method of approximation used in this analysis. 

Figure 7(b) shows results for a slightly stronger field which provides 
more radial gain at the cost of greater vertical defocussing. 


$7. THE ELECTRONIC COMPUTER AT ROTHAMSTED 


To check the approximate calculations a numerical solution of eqns. (7) 
and (8) was carried out on the electronic computer housed in the 
Statistical Department at Rothamsted. This machine is the ‘401’ 
prototype model built by Elliott Brothers under contract from the 
National Research Development Corporation (Lipton 1955). 

The 401 is a two address serial machine with a word length of 32 
binary digits, one of which is a sign digit—this corresponds to about 
95 decimal digits. The main store is a magnetic disc of 23 tracks, 
each track holding 128 words. All the usual operations and functions 
are available with the exception of division for which a sub-routine is 
used. Information is fed into the machine by 5-hole punched tape. 
Output is by electric typewriter or teleprinter punch. 


§ 8. MerHop oF SoLuTION 


A modified Runge-Kutta process was used (Wilkes, Wheeler and 
Gill 1951). Equations (7) and (8) were written as four simultaneous 
linear equations and to avoid overflow p, z, 2 were each scaled down by 
a factor of 10 and p by 100. (In the machine a number, ¢, lies in the 
range (—l<t<1). Thus the equations to be solved were 


u=10e u=10x 
v=y v=y 
. 7 Rod p<0, ‘ >0 
x=—0-0955u x= —0-955u(0-1—0-152w)—0-13v2 
y= —0-045v y= — (0-045-+2-6u)v j 
where 


10u=p, 10v=2z, l0y=dz/d0, 100r=dp/d@. . . . (24) 
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The scaled equations for the new values of w and v after regeneration 
(assuming p>0) were 


Wout = Uin— L'u—10Vu?+ 100Wu8 


Vout= Vin + Lv-+ 20 Vuv— 300 Wwv eee (25) 
for a field given by 
r 
O, (7) Tet Woe Wie. ee ee (26) 
0 


The step length of integration was 4°. Results were printed out at 
values of 6 of 328° and 360° during each revolution (the values at 360° 
were before regeneration). The information printed was the values of 
p, 2, dz/d0, dp/d0, the maximum value of z reached go far and g, the value 
of 6 at which p becomes negative, as in fig. 2. 


§ 9. THE PROGRAMME 


In order to apply the Runge-Kutta process three main routines are 
required. One of these is the basic integration routine which in turn 
uses an ‘ auxiliary’ routine four times during each step. The purpose 
of the auxiliary is to calculate w, v, x, and y from the values of the variables 
previously computed by the basic routine. Four such operations are 
required to advance the integration one step. In addition a ‘ master’ 
routine coordinates the computations, keeps a count of the number of 
steps and prints the results when necessary. 

A complication was introduced by the action of the regenerator. 
This difficulty was overcome by, in effect, restarting the computation 
at the end of each cycle of 90 steps: the values of p, z at this stage and 
the values of dp/d0, dz/d@ were used as new initial conditions. The 
master routine ensured that this operation occurred at the right value of 6, 
and took place only if p>0. 

Originally the programme was written so that the relevant defining 
eqns. (7) or (8) became operative as soon as p changed sign. This 
procedure was eventually altered and the range of negative p was crossed. 
in one step in the following way. 

If 6, is the value of 6 at which p becomes negative then it is permissible 


to jump immediately to 6,,-, where 
Oy = pee e010: 9551? —0,,4- 184" 
with 
Ui = — Uns Uy = — Un 
vy, =0-7859v, +2:9160,,, 
y= —0-13110, +0-78590,,. 
This modification considerably reduces the time taken to complete a 


whole revolution. 
The Runge-Kutta has a truncation error in one step of order h°, where 


h is the step length, but rounding off errors accumulate at a rate inversely 
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proportional to the number of binary digits held. The results are 
certainly correct to 3 decimal places. 

Integration round a complete revolution with the modified programme 
described above took about 1 minute, including two rows of printing, 
and about 12 revolutions were obtained for each set of initial conditions. 
The solutions were computed over a period of four months, the work 
being carried out only during short intervals when the machine was not 
being used for its more usual statistical computations. 


§ 10. REesuLTS 
Several starting values of the variables were used. One complete set 
of results served to indicate the effect for several different sets of initial 
condition, for each set of values printed at 360° in every cycle could be 
regarded as a new set of initial conditions. Thus we have the tracks of 


Fig. 8 


Regenerator field shapes H q for cases A, B, C, D, E. 


particles starting with initial values p0-5 in, po, z2=0-5 sin o. 
ia 4 ¥ = = . . < i . . ‘ me 
2=0 5w cos « for any phase «; particles with smaller initial vertical 
displacements have smaller subsequent displacements approximately in 
proportion. : 
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Table 1 lists the regenerator field shapes (26) which were tested, and. 
for easy comparison these are shown in fig. 8 and their performance in 
figs. 9 and 10. 


Table 1 
Field cp 
A 0-2 
B 0-2 
Cc 0-2 
D 0-2 
E 0-2 


Fig. 9 


The increase, p,+1—Py» Of successive radial displacements at 32° before the 
regenerator. 


Figure 9 shows the radial separation p,-1—Pn of successive passages 
past azimuth 328°, that is 32° before the regenerator, together with the 
phase angle g of the motion. Figure 10 shows the maximum vertical 
displacement experienced anywhere by any particle with an initial 
amplitude of } in., that is the height of the beam. 


1278 K. J. Le Couteur and 8. Lipton on Non-Linear 


The vertical displacements at 328°, where an outlet channel might be 
placed, are usefully smaller than the maxima shown in fig. 10 because 
where p is large and positive the vertical focussing forces are very strong. 

Clearly curves A and B give a satisfactory combination of limited 
vertical expansion and adequate radial separation, im reasonable 
quantitative agreement with the approximate results for these fields 
shown in fig. 7(a), (b). Fields C, D, E were used to study the effect of 
deviations from the shape indicated theoretically and are less satisfactory 
than A or B for they lead to dispersal of the beam when p exceeds 3 in. ; 


Fig. 10 


fo) 1 2 3 4 5 6 PrP 7 in, 8 


Maximum vertical displacement for any particle with less than 0-5 in. initial 
amplitude of vertical oscillation. 


some particles go off vertically and retreat inwards while others continue 
outwards. For C at first sight this is surprising since it les between 
A and B. The reason seems to be that for C, D, E the gradients dH ,/dr 
are less than for A and B and for p>3 in. are too small to satisfy the 
vertical stability criteria established in § 5. 

It is noteworthy that fields A and B produce radial displacements 
greater than 7 in., which carry particles beyond the point of free escape 
from the cyclotron. 


§ 11. CONCLUSION 


The theoretical analysis is easily adapted to any synchrocyclotron 
and shows that a non-linear regenerative deflector system can extract 
the beam with a higher energy than that now working at Liverpool. 
To get an intense external beam the particles must finally leave the 
machine through a magnetic channel, but since the deflector carries 
particles out to large radii the channel need not produce much field 
reduction. The resulting possibility of using a short magnetic channel 
of large aperture should lead to still greater efficiency than has already 
been achieved. 

Whereas the linear apparatus requires two deflector fields, known as 
peeler and regenerator, the non-linear apparatus requires only one with 
a curvature in the sense which arises naturally when the field is produced 
by blocks of ions fixed in the synchrocyclotron. 
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The mathematics has some interest apart from the application to 
beam extraction, for the difficulty of handling the coupled non-linear 
equations (7) and (8) in the presence of similar periodic perturbations 
is at present one of the main obstacles to progress in the design of new 
particle accelerators of the fixed field type. The reasonable agreement 
found between results of a simple approximation and of extensive 
numerical computation suggests that these equations may be more 
tractable than is generally assumed. 

These proposals were primarily intended for the Liverpool and 
C.E.R.N. synchrocyclotrons, but are now being applied to several others. 
The expected results have already been obtained at Chicago, by Dr. Crewe, 
using the deflector fields shown in fig. 7(a). 


Ace EN eT xX 


As the precise characteristics of the deflected beam are important, 
table 2 gives particle coordinates at the last few successive passages 
past azimuth 328° for field A and various initial conditions. p vanishes 
at azimuth g+2° and then at azimuth g+184+2°, which is the virtual 
source of particles swinging outwards. 


Table 2 
Revolution p(in.) 2(in.) dz/d0 dp|d0-—“2nae g(-) 
initial 0-5 0-5 0 0 
10 2-57 0-53 —0-40 0:64 0-86 58 
if 3-57 =0-77 0-62 1-11 —1-02 58 
12 5-80 0-82 Bae 2.48 158 58 
13 +10 
initial 0-5 0 0-125 0 
10 2-17 0:39 —0-07 0-51 0-72 58 
11 2.96 0-18 —0-24 0-90 0-72 62 
12 4-4] —0-55 0-67 1-69 0-86 58 
13 8-40 0-82 ==[:55 4-71 235° 54 
initial 0-5 0-2 0-1 0 
10 2-09 0-76 —0-30 0-38 0-86 58 
11 2-89 —0-26 —0-00 0-79 0-86 58 
12 4-36 —0-28 0-46 1-66 0:86 58 
13 8-50 0-63 Apa 4-94 1:82 54 
ee 0-5 0 0-0125 0 . 
or I 2-14 0-054 —0-020 0-63 0:07 62 
12 2-88 OL —0-006 0-96 0-07 62 
13 4-33 —0-032 0-043 1:75 0-07 62 
14 8-94 0-067 —0-079 5-62 0:16 54 
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SUMMARY 


Measurements have been made on the single (large angle) and multiple 
coulomb scattering of electrons and positrons in nuclear emulsions. The 
results obtained for single scattering are in general agreement with the 
theoretical expectations, and appear to confirm the prediction of a 
difference between the nuclear scattering of electrons and positrons at 
10 Mev, although the data are insufficient to permit detailed comparison. 
The results for the multiple scattering show an unexplained difference 
between the two types of particle. The number of positrons observed to 
annihilate in flight is in agreement with the theoretical predictions of 
Dirac (1930). 


§ 1. SINGLE SCATTERING 
(i) Introduction 


THE single or large angle scattering of relativistic electrons has been 
studied by several workers at energies less than 4 Mev, by Lyman, Hanson 
and Scott (1951) at 15-7 Mev, and by Hofstadter et al. (1954) at higher 
energies. In contrast with the widely varying results of earlier experi- 
ments the more recent work at the lower energies is in general agreement 
with theoretical predictions. Apart from the work of Fowler and Oppen- 
heimer (1938) at about 11 Mev, in which only a few deflections of over 14° 
were observed, the single scattering of positrons has been studied only 
at low energies (<2 Mev ; Howatson and Atkinson 1951). The nuclear 
scattering of electrons and positrons with energies in the range from 
4-15 Mev provides a test of the validity of the relativistic scattering 
formula for point charges, but can give little information about the 
nuclear charge distribution. In the present experiment the single 
scattering cross sections at 10 Mev for electrons and for positrons were 
examined under conditions which enabled a direct comparison to be 
made between them. 


(ii) Haposure of the Emulsions 


The Glasgow H.T. set was used to provide the 10 Mev positrons and 
electrons by means of the nuclear reactions *Li(p, y)*Be and “Li(d, p)*Li. 


* Now at Ewell County Technical College, Surrey. 
+ Communicated by Professor P. I. Dee, F.R.S, 
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In the first case both types of particle were produced from the 17:6 Mev 
gamma-rays by pair production in lead, and in the second case a large 
flux of electrons with a maximum energy of about 12 Mev was obtained 
from the decay of the short-lived *Li nuclei. All the positrons were 
produced from the first reaction, and electrons were obtained by both 
methods. The greater yield from the 7Li(d, p)*Li reaction allowed much 
shorter times of exposure than were possible in the case of pair production. 
The particles produced in these reactions were analysed by a double 
focusing magnetic spectrometer which had a resolution rather better than 
5%, for the absolute value of the energy, and were then allowed to pass 
into 400 » Ilford G5 emulsions. The magnetic field was measured at 
intervals during the exposure of each plate and the corresponding energies 
in all exposures were close to 10 Mev. ‘The final results have all been 
normalized to this energy. 
(iii) Analysis 

The measurements were made on Cooke M. 4000 microscopes using 
x15 eyepieces and x45 oil immersion objectives. The plane of the 
emulsion was inclined at ~9° to the beam of particles which extended 
over the central inch of each three-inch plate. Only tracks from the 
central portion of the beam were used for measurement since a check, by 
multiple scattering measurements, on the energies of particles distributed 
over the beam revealed a few rather low-energy particles on the fringe of 
the beam in one exposure. Apart from this limitation the acceptance 
criteria for a track were 


(a) that it should lie within 10° of the mean beam direction, 


(6) that it should have a minimum length of 500, for surface to 
surface tracks and of 700 for surface to glass tracks. 


Suitable tracks were followed from their point of entry into the emulsion 
to a maximum length of 2000. For all scattering angles having a 
horizontal projection greater than 10° measurements were made of the 
horizontal projection of the angle of scattering and of the dip angle of 
the primary and scattered tracks. From these measurements the spatial 
angle of scattering was calculated in each case. Track lengths were 
measured using an eyepiece scale and corrected for the effect of dip. 
The experimental results were corrected for effects due to 


(a) the rejection of scattering angles with horizontal projections <10°, 

(b) the rejection of tracks having lengths less than 500 or 700 microns. 
(iv) Results 

380 cm of electron track and 520 cm of positron track were scanned 

giving total numbers of deflections > 25° of 136 for electrons and 163 for 

positrons. The resulting probabilities for scattering in the emulsion, 


after the application of the correction factors mentioned above, are 
shown in table 1. The errors indicated are standard deviations, 
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Table 1 also gives the anticipated theoretical values for the probability 
of scattering as obtained from the calculations of Feshbach and others 
(Feshbach 1952, McKinley and Feshbach 1948, Acheson 1951) for electrons 
and Yadav (1952) for positrons. In obtaining the theoretical values we 
have evaluated the sum 2;N ,o,(0,—9.) where o,(0,—6,) is the scattering 
cross section of the ith element for the angular interval (6;—0,), N; is 
the number of atoms of the ith element present per cm® of the emulsion, 
and where the sum extends over all the elements present in the emulsion. 
Since oc; is proportional to z,2, where z; is the atomic number of the ith 
element, nearly all the scattering is due to the silver and bromine content 
of the emulsion. The effect of finite nuclear size has been included in the 
calculation of the cross sections by using the results of Acheson (1951) 
and of Elton and Parker (1953), but the effect is small (~7%), even for 
the highest angular interval. 


Table 1. Probability of Scattering per cm of Track Length 
in the Emulsion 


Angular interval 25°-35° 35°-50° 50°-180° 
Measured. 0-275-40-031 | 0-120+0-020 | 0-110 +0:019 
Electrons =_|_—————_|-— 
Theoretical 0-280 0-150 0-102 
Measured 0-220-+0-024 | 0-115+0-017 0-066 0-013 
Dien | aa Saray 


Theoretical 0-203 0-095 0-054 
pw Ee ee 


The statistical errors associated with the experimental data are too 
large to allow a detailed examination to be made of the differences in the 
scattering probabilities for electrons and positrons as a function of 
scattering angle. The results do, however, appear to confirm the 
theoretical prediction of a difference between the nuclear scattering of 
electrons and positrons at 10 Mev due to the differences in the spin—orbital 
forces on the two types of particles, and also show that the theories of 
single coulomb scattering at this energy predict cross sections of the right 
magnitude. 


§ 2. MULTIPLE SCATTERING 


The multiple coulomb scattering of electrons and positrons in nuclear 
emulsions has been studied by several workers (Corson 1951, Gottstein 
et al. 1951, Bosley and Muirhead 1952, Husain 1955, Williams, private 
communication). 

Direct comparisons between the multiple scattering of electrons and 
positrons with energies of about 2 Mev in argon have been made by 
Groetzinger, Humphrey and Ribe (1952), and by Cusack and Stott (1955), 
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using a cloud chamber. A study of the scattering in emulsions of both 
types of particle under the same conditions has been made by Corson 
using particles of 196 Mev and by Williams using particles of 1, 15 and 
2-5 mev. The present work enables a direct comparison to be made 
between the scattering of electrons and positrons with an energy of 10 Mev. 

Earlier work in this laboratory with a magnet of poorer resolving power 
had indicated a difference between the multiple scattering of electrons and 
positrons. .In the present work tracks were measured on the same plates 
as were used for the single scattering experiment, and the relative 
energies of the two types of particle should be correct to about 2%. 

Scattering measurements using the coordinate technique developed by 
Fowler (1950) were made on about 9 em of electron track and about 
10 cm of positron track, and from these results values of the scattering 
constant K for the emulsion were calculated from 


Keates (>) ay 


cHorp 2 \ ¢ 


where <¢)cHorp is the mean angle of scattering measured by the 
coordinate method, p, Be are the momentum and velocity of the 
scattered particle, z is the ratio of its charge to that of the electron, 
and ¢ is the cell size in microns. K is expressed in the usual units of 
degrees < Mev X (100 p)~¥/2. 

The measured values of AK were corrected for the effects of spurious 
scattering due to irregularities in the microscope movement and to 
observer errors (‘ noise ’), for the inclination of the tracks to the coordinate 
axes used for measurement, and for distortion of the emulsion. Except 
for the very small cell sizes where the effect of ‘ noise’ is appreciable 
(‘ noise’ correction for 25 cell is ~10°%), all these corrections are 
small. 

The corrected values of K are plotted in fig. 1(a) as a function of cell 
size. A value of the scattering constant less liable to be affected by 
variations in the small number of large deflections in the distribution is 
obtained by omitting angles greater than four times the median value. 
These ‘ cut’ scattering constants K, are plotted in fig. 1(d). 

In table 2 the results obtained from the present work are compared 
with the theoretical predictions and with results obtained by other 
workers. The theoretical results quoted here are those of Moliére (1948), 
and take no account of the spin of the particles. The magnitude of the 
difference in scattering to be expected from this effect has been calculated 
by Mohr and Tassie (1954) ; their results predict a difference in scattering 
between positrons and electrons of less than 1°, under the conditions of 
the present experiment. The results quoted in the table do not include 
those of Williams (loc. cit.) for the scattering of 1, 1-5 and 2-5 Mev 
electrons and positrons in nuclear emulsions, which are presented in 
a forthcoming paper, where a difference is also indicated. The results 
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from the cloud chamber measurements are not included in the table, 
since they do not yield values of a scattering constant directly comparable 
with those obtained using nuclear emulsions. The measurements of 
Groetzinger, Humphrey and Ribe (loc. cit.), at 1-5—2-5 Mev show a difference 
of about 10% between the scattering of electrons and positrons compared 
with a theoretical prediction ~1%. Cusack and Stott (loc. cit.), at 


0-4 Mev find no significant difference between the scattering of the two 
types of particle. 


Fig. 1 


35 
ELECTRONS © oH] ELECTRONS ae 
POSITRONS * POSITRONS 
THEORETICAL —— THEORETICAL ——— 


ie) 100 200 300 400 Oo 100 200 300 400 
CELL SIZE IN MICRONS CELL SIZE IN MICRONS 


(a) (6) 
Table 2. Values of the Scattering Constant K for 100 p Cell Size 


Author Positrons Electrons 


Present work, 10 Mev 24-9+0:5 28-6+0°5 


Gottstein et al., 105 Mev 25:-8+0°6 


Corson, 196 Mev 25:0 40:8 


Bosley and Muirhead, 9-19 Mev 


Husain, 6-6 and 11-5 Mev 


Theoretical (Moliere) 


§ 3. ANNIHILATION OF PosITRONS IN FLIGHT 


The positron plates showed 18 examples of particles which appeared 
to annihilate in flight. The number calculated from the theory of Dirac 
(1930) is 22. 


oe 6 Ae 
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ABSTRACT 


141Ce has been aligned by the method of magnetic hfs (alignment) 
in a crystal of cerium ethyl sulphate. From measurements of the 
angular distribution and the plane of polarization of the y-radiation it 
is concluded that the y-ray transition is principally M1 with a small 
admixture of E2 (5=+0-08-L0-02). Values of 0-75+-0-20 and 0:66-L0-16 
nuclear magnetons are found for the nuclear moment of 141Ce, depending 
on whether the f-transition involves a spin change of | or 0. 


§ 1. INTRODUCTION 


THe method of nuclear orientation applied to the investigation of nuclear 
properties by observing the directional and polarization effects of the 
y-radiation from oriented radioactive nuclei has been described in earlier 
papers (Bleaney et al. 1954, Bishop et al. 1955 a). It was indicated in these 
that for a given cascade of B- and y-ray transitions and for a given assign- 
ment of spins to the associated levels the anisotropy of a particular 
y-ray depends on the multipole order of the transition, on the character 
of previous transitions and on the original degree of orientation. It was 
also pointed out that quite small coherent admixtures of other multipoles 
in the y-ray transition can produce an angular distribution radically 
different from that of the pure transition. Atthesame time this admixture 
will alter the degree of plane polarization of the radiation but to a different 
extent. For instance, it is possible to destroy the y-ray anisotropy 
completely and yet leave a large polarization. Conversely it is possible 
from simultaneous measurements of the y-ray anisotropy and polarization 
to determine the mixing ratio in mixed y-ray transitions. This method 
of studying mixed transitions has now been applied to the isotope 1*Ce 
which was oriented using the method of magnetic hfs (alignment) ; this 
isotope was first aligned by Ambler, Hudson and Temmer (1955). The 
present experiment differs from that of Ambler ef al. in two respects. 
Firstly a crystal was used which gave nuclear alignment along an axis in 
contrast to alignment in a plane. Secondly the polarimeter was employed 
which showed that the y-ray contained a mixture of multipole radiations. 
ie meee ee SS ee 


* Communicated by the Authors. 
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§ 2. Decay ScHEME OF 141Ce 


The principal features of the decay scheme of ''Ce are well known and 
are shown in fig. 1 (Jones and Jenson 1955). The spin of the ground 
state of “Pr has been measured as 5/2 and from its nuclear moment this 
is identified as a d 5/2 state on the nuclear shell model. The assignment 
of f 7/2 to the ground state of 11Ce is made on the basis of comparative 
half-life values for the f-transitions which suggest that both are first 
forbidden with i;=0or +1. Internal conversion (Heath 1952) and lifetime 
measurements (de Waard and Gerholm 1955) on this y-ray indicate that 
it is principally M1 and this appears to be confirmed by the nuclear 
orientation measurements of Ambler et al. These results can be most easily 
understood by an assignment of g 7/2 to the first excited state of “1Pr 
(see fig. 1). 


145 KEV 


The decay scheme of Cerium 141. 


Temmer and Heydenburg (1954) have shown that the 145 kev state is 
weakly excited by Coulomb excitation of 1Pr and we conclude that at 
least part of the y-ray transition must be electric (presumably E2). By 
comparing their figures on the intensity of this radiation with that obtained 
in '*"T'a we can estimate the partial half-life for E2 decay as approximately 
3x 10-8 sec. Direct measurement of the half-life of this state for y-ray 
emission by de Waard and Gerholm (1955) gives a value of 2-4+-0-3 x 10-9 
see and we therefore conclude that the ratio of the intensities K2/M1 is 
about 10-8. Since the half-life estimate from Coulomb excitation is very 
approximate this value for the ratio of the intensities is only correct to 
within an order of magnitude. 
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§ 3. Previous NUCLEAR ORIENTATION RESULTS 

The experiments of Ambler et al. were made using a crystal of cerium 
magnesium nitrate. In this crystal the cerium nuclei become oriented in 
such a way that the nuclear magnetic state with 1,=-+} is preferentially 
populated, i.e. the nuclei tend to align with their spins in a plane perpen- 
dicular to the z-axis. From the observed y-ray angular distribution and 
its dependence on temperature these authors conclude that the transition 
is M1 and also that the spin of the 145 kev state is 7/2*. To determine 
the nuclear moment of 14!Ce it is necessary to know the hfs field for cerium 
in this crystal. The hfs for cerium has not been measured experimentally 
since no stable isotope with a moment different from zero exists. However, 
using calculated values for the hfs separation they deduced that the nuclear 
moment is 0-16 nuclear magnetons. This is smaller than that to be 
expected for a nucleus whose moment arises from a single f 7/2 neutron ; 
other nuclei with f 7/2 neutron configurations have moments of about 
1 nuclear magneton (Walchli 1953). This small value could arise from 
two possible causes. Firstly, there might be a reduction in the hfs field 
due to interaction between the cerium ions, not taken mto account in the 
calculation. This could give rise to a smaller degree of orientation of 
141C¢e@ and hence lead to a smaller value of nuclear moment. This is 
discussed by Ambler e¢ al. Secondly, there might be a reduction in the 
y-ray anisotropy either due to a perturbation of the first excited state of 
141Py through interaction with extranuclear fields, or to the presence of 
a coherent admixture of E2 radiation in the Ml transition. It seems 
reasonable to suppose that in planar alignment, where the lower levels are 
relatively closer together, the nuclei would be more likely to be perturbed 
by extra-nuclear fields than in axial alignment. 


§ 4. Dusign or EXPERIMENT 


‘To avoid those difficulties which may be present in the cerium magnesium 
nitrate crystal we decided to align this nucleus in a crystal where axial 
alignment should be obtainable. This is a more favourable case for 
obtaining orientation as the splitting between the lowest levels is in general 
larger (Bleaney 1951). Because of lack of knowledge of the hyperfine 
structure field in cerium one must be guided by g-values in choosing the 
crystal. In cerium ethyl sulphate g,=3°9, g,=0-2 (Bogle, Cooke and 
Whitley 1951) which implies that A is greater than B and hence alignment 
along the axis would. be obtained. However one can only reach tempera- 
tures of about 0-08°K (Cooke, Whitley and Wolf 1955) and consequently 
only low degrees of orientation are obtained. If the cerium is diluted 
with lanthanum it should be possible to reach much lower temperatures 
but unfortunately in the dilute salt g,, is less than g, and hence the lowest 
nuclear state has I1z—=-+}. Thus again alignment in a plane would be 
obtained. 

The difficulty was overcome by replacing all but approximately 1% of 
the cerium which included the radioactive isotope 1*41Ce by neodymium. 
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In an auxiliary experiment we had shown that concentrated neodymium 
ethyl sulphate can be cooled to about 0-024°K by adiabatic demagnetization 
and this salt has already been used to align 147Nd nuclei (Bishop et al. 
1955 b). The hfs of cerium in this crystal is not expected to be markedly 
different from that in concentrated cerium ethyl sulphate. Expressions 
for A and B in the spin Hamiltonian of Abragam and Pryce (1951) may be 
obtained from the treatment of Elliott and Stevens (1951, 1952). The 
value of (1/r3) where r is the radial coordinate of the 4f electron with 
respect to the nucleus has been given by Bleaney (1955) as 32-5A7%. 
In order to evaluate A and B explicitly knowledge of the Jz value of the 
lowest state is necessary. As will be shown later, axial alignment of 
cerium in this salt is obtained and therefore the state with Jz;=--5/2 lies 
lowest. The corresponding values of A and B are 


A=0°108 (uy/Z) cem-}, 
B=0-:002 (y/Z) em}. 


A pair of these neodymium ethyl sulphate crystals weighing about 
2 g was mounted in the low temperature cryostat with their crystalline 
axes horizontal. These contained about 13 microcuries of 144Ce, prepared 
by irradiation of stable cerium oxide in the Harwell pile. This arrange- 
ment is convenient for the y-ray distribution measurement since it allows 
the counters to be mounted in the horizontal plane. For the neodymium 
ions 


Qu=3'61, g,=2-05 (Bleaney, Scovil and Trenam 1952). 


Since these g-factors are of the same order of magnitude the direction of 
the large magnetic field used in cooling by adiabatic demagnetization is 
relatively unimportant. The magnetizing field was applied vertically 
along the g, direction. Using values of H/T of 20 kilogauss/°K magnetic 
temperatures 7’'® of about 0-024°K were reached. For these crystals it 
was found that the temperature remained below 0-:03°K for about 20 
minutes. 


§ 5. ANGULAR DISTRIBUTION AND POLARIZATION OF THE y-RADIATION 


On reaching the low temperature the magnet was removed and a pair 
of Nal (Tl) scintillation counters was mounted round the cryostat. 
Pulse height selection was used to discriminate against x-radiation and 
radiation degraded by scattering in the apparatus. The angular distri- 
bution was observed at a series of temperatures reached by demagnetizing 
from different values of the magnetic field. In fig. 2 is shown the depen- 
dence on temperature of the y-ray anisotropy «; «¢ is defined ‘as 
1—W(0)/W(z/2) where W(6) is the intensity of the y-radiation emitted at 
an angle @ with the axis of orientation. It is seen that the anisotropy is 
negative, whereas Ambler et al. found positive values. We therefore 
conclude that alignment along the crystalline axis is obtained as expected. 
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Measurements were also made on the plane polarization of the y-radia- 
tion emitted in the plane perpendicular to the alignment axis. These 
were made with the polarimeter which has been described elsewhere 
(Bishop and Perez y Jorba 1955) and makes use of the dependence of the 
differential Compton scattering cross section on the plane of polarization 
of the y-radiation. Polarimeter measurements were confined to the lowest 
temperature region where the effect was largest and these showed the 
correct general variation with temperature. Defining P’ as the ratio of 
the numbers of quanta scattered in the directions perpendicular to (V,) 


Fig. 2 
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The dependence of y-ray anisotropy (c) from aligned Cerium 141 nuclei on Lie. 


and parallel to (N,,) the axis of alignment, we find that for 1/7 © of 40 P’is 
1:58. This ratio P’ is connected with the polarization of the radiation ‘Pe 
by the relation 

P=(1—P’\(R4ED/U+P')\(R—1). 


P is defined as (I, —1,)/(Z, +41) where J, and /,, are the intensities of the 
radiations polarized with their electric vectors perpendicular to and parallel 
to the axis of alignment, and R is the analysing ratio for the polarimeter 
and depends on its geometry and the y-ray energy. In this case R is 
believed to be about 6. This value has been estimated from the known 
behaviour of the polarimeter. Experiments to test the validity of this 
assumption are at present in progress. 
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§ 6. INTERPRETATION OF RESULTS 


The intensity of radiation W(0, 4) emitted at an angle 6 to the axis of 
orientation from a system of nuclei can be derived from the expressions 
given by Biedenharn and Rose (1953) for polarization—direction correlation 
between successive y-rays. ¢ is the inclination of the electric vector to the 
plane containing the axis of orientation and the direction of propagation. 
The parameter F',(L,j,,j) which describes the first y-ray is replaced by the 
orientation parameter B,U,; B, describes the orientation of the initial 
state (in this case the ground state of 4Ce) and U,, depends upon the spins 
involved in the unobserved transitions preceding the y-ray (in this case the 
B-transition)(Gray and Satchler 1955). For low degrees of orientation 
where B, is the only significant term and for dipole transitions with 
admixtures of quadrupole radiation this expression reduces to 

W(6, d)=1+ BU oF Q(1, Lp, 11) P2 (cos 0)(1+36 + 386?) 
+4o B,U,F,(1, Io, 11) P.? (cos @) cos 26 (1—«6+ 3882). 

P.(cos 6), P.2(cos 6) are Legendre and associated Legendre polynomials. 
5 is the ratio of the amplitudes of quadrupole to dipole radiation. It 
should be noted that ¢ as defined above is different from that used by 
Biedenharn and Rose, being complementary to their angle. 

a= 2(—)b-P-1.4/15 (21,2 D Ge (1, 1p, 1q)/3Fs gy) 

Beeks (2; Lada ion od eee eee 
F, and G, are functions tabulated by Biedenharn and Rose (1953). 
o takes the values +1, —1 for magnetic and electric dipole transitions 
respectively. 

Using this expression one can first of all deduce from the signs of P and 

« that the transition is principally magnetic (either M1 with a small 
admixture of E2 or M2 with a small admixture of El). The second 
choice is unlikely both because the internal conversion coefficient and the 
lifetime of the state favours an M1 assignment and also because the value 
of 6 required leads to rather small predicted values for the anisotropy at 
any temperature. 

From the value of +1-58 for P’ at 1/7’© =40 we find that P is —0-30. 
Combining this with the value of —0-157 for e at the same temperature 
we find that 6 is +-0-08+-0-02 for the E2/M1 mixing ratio. This result 
depends only on the spin assignments 7/2 and 5/2 for the initial and final 
states of the y-ray transition and is independent of any possible perturba- 
tion of the states. The experimental error arises from the statistical errors 
in P’ and ¢ and in uncertainty in the value of R. The theoretical curve 
shown in fig. 2 has been calculated using this value of 6 and has been 
adjusted to pass through the point e=—0-066 at 1/7® =25. This point 
was chosen because for lower temperatures the difference between 7’ ® 
and 7’ although unknown may be considerable and for higher temperatures 
the anisotropy becomes too small to be measured accurately. The devia- 
tion of the points from the curve in the low temperature region is attributed 
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to differences between 7’© and 7’. With this value for 5 and a value of 
0-066 for at 1/7’ © =25 we find that B,U,F, (1, Ip, 1,) is 0-0810-017. In 
this way we find that B, is 0-302-+.0-066 if the £-transition involves a 
spin change 7g—1 and 0-246--0-054 if the spin change 7g=0. From the 
corresponding values of 8(=A/2k7') we find that A is 0-0231+-0-0059 cm~? 
(ig=1) or 0-0201+0-0049 cm™? (tg=0). From the calculated values 
for AI/w (see § 4) we find that the nuclear moment of '1Ce is 0-75--0-20 
(tg=1) or 0-66+0-16 (1g=0) nuclear magnetons. The large errors here 
arise principally from uncertainty in the value of & and include an allowance 
of 10% for the possible difference between absolute and magnetic 
temperatures. - 

These figures for the nuclear moment are markedly higher than that 
obtained by Ambler et al. Part of this difference can be attributed to the 
effect of E2/M1 admixture whose neglect would lead these authors to 
conclude a smaller value for this moment. If one takes our figure for 6 
and recalculates their value for the nuclear moment a figure of 0-27 
nuclear magnetons is obtained, which is still considerably smaller than 
our value. 

The remaining difference can arise from two possible effects. Hither the 
calculation for the hfs in cerium magnesium nitrate is in error due to 
interaction between the cerium ions which has not been taken into 
account or the 145 kev state in !!Pr experiences a perturbation during its 
lifetime due to interaction with extra-nuclear fields. It seems plausible 
that such a perturbation could have a greater effect on nuclei oriented in a 
plane, as in the double nitrate salt where the lower hfs levels are separated 
by about 0-002 cm}, than on nuclei oriented along an axis as in the ethyl 
sulphate where the separation is 0-01 cm~!. | In this case we can estimate 
the magnitude of the perturbation if it is assumed to be negligible for the 
ethyl sulphate. If we define an attenuation coefficient A by 

W(0)=1+A B,U2F 2 (1+3a5) P2 (cos 4) 
where A takes the value 1 for no perturbation and 0 for complete destruc- 
tion of the orientation, then it is found for the double nitrate that A is about 
0:5. In view of the uncertainty as to whether a perturbation is also present 


in cerium ethy! sulphate we must conclude that our values for the nuclear 
moment represent lower limits. 


§ 7. CONCLUSIONS 


It is concluded from these experiments on the orientation of 41Ce in a 
crystal of cerium ethyl sulphate that the 145 kev y-ray transition excited 
in its decay is predominantly M1 with a small E2 admixture as is suggested 
by the Coulomb excitation experiments (see § 2). The value of the mixing 
ratio 5 is 0:08-0:02. The nuclear moment of 41Ce is 0:75-0:20 or 
0-66-40-16 nuclear magnetons depending on whether the f-transition 
involves a spin change of 1 or 0 and assuming that there is no perturbation 
of the 145 kev state. There is a difference between the nuclear moment 
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obtained in these measurements and in those of Ambler e¢ al. Part of 
this can be accounted for by the presence of the E2/M1 mixture and the 
remainder must be attributed either to a perturbation of the nuclear state 
in the cerium magnesium nitrate, or to the existence of some other form of 
coupling not taken account of in the calculation of B in the double nitrate 
crystal, possibly due to the presence of neighbouring cerium ions. 
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SUMMARY 

Accurate measurements have been made of the lattice spacings of 
dilute solid solutions of titanium, vanadium, chromium, manganese, 
cobalt and nickel in «-iron. In all cases the formation of the solid solution 
results in an expansion of the iron lattice, and for equal atomic percentages 
of the solutes, the expansions are in the order Ti>V>Cr and Co<Ni<Cu. 
Manganese behaves anomalously, its effect being nearly the same as that 
of chromium. With titanium and vanadium as solutes there are negative 
deviations from Vegard’s law, whilst for the remaining solutes the 
deviations are positive. These results, together with previous work on 
nickel alloys, are discussed, and some of the effects are ascribed to 
exchange repulsion between nearly filled d shells. 


§ 1. INTRODUCTION 


Tue object of the present work was to measure the lattice spacings of the 
solid solutions in «-iron formed by the elements from titanium to nickel, 
and to compare the results with those obtained by Hellawell and Hume- 
Rothery (1954) for the solid solutions in ruthenium formed by the corre- 
sponding elements in the Second Long Period, and with those obtained 
by different workers for solid solutions in nickel. 


§ 2. EXPERIMENTAL 
(a) Metals Used 

The iron used was special high purity H iron prepared at the National 
Physical Laboratory, and obtained from the Pure Metals Committee of 
the Department of Scientific and Industrial Research ; we have to thank 
Dr. H. Sutton for his help in this connection. The analysis indicated a 
purity of the order 99-969%, the chief impurities being (in weight per- 
centages) carbon 0-006, manganese <0-005, nickel 0-006, copper 0-004, 
and aluminium 0-004. The titanium used was I.C.I. Grade A iodide metal 
whose purity was of the order 99-7%, the principal impurity being iron and 
therefore harmless. The vanadium used was supplied by Messrs. Johnson 
Matthey & Co. Ltd., and was in the form of beads of purity 99-7-99-8%,. 
For the work on iron—chromium, iron—manganese and iron—nickel alloys, 
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use was made of bars of specially pure alloys prepared for work at the 
National Physical Laboratory, and we must thank Dr. N. P. Allen for 
presenting these, some of which were used as provided, whilst others were 
remelted with the high-purity iron. For the cobalt alloys, use was made 
of a master alloy prepared by arc-melting by Messrs. Johnson Matthey & 
Co. Ltd., and the purity was of the order 99-93%, the chief impurities 
being nitrogen, oxygen and nickel. 


(b) Preparation of Alloys 


The iron-titanium alloys were prepared in an argon arc-melting furnace, 
the earlier melts being carried out at the R.A.E. Farnborough, by the 
kindness of Mr. L. G. Carpenter. Later alloys were prepared by arc- 
melting in our own laboratory. Some of the iron—vanadium alloys were 
prepared by arc-melting at the R.A.E., whilst all other alloys were 
prepared by melting in an H.F. induction furnace under low pressure 
of hydrogen, with full precautions of evacuation and degassing. The 
crucibles used were of recrystallized alumina or thoria, and the alloys 
were allowed to solidify in the furnace. For the cobalt alloys the 
hydrogen was pumped off after melting, and solidification took place 
im vacuo. 

Specimens from the ingots were homogenized by heating for 3-5 days 
~ at 1200°-1250°c in alumina collars contained in silica tubes which were 
either evacuated or filled with argon to an extent producing a pressure of 
about 1 atmosphere at the annealing temperature. After homogenization 
the alloys were air-cooled, except that the iron—nickel and iron—manganese 
alloys were quenched in order to prevent decomposition on cooling. The 
N.P.L. alloys had already been homogenized, but were re-annealed and 
quenched, except for the iron—chromium alloys which were used as 
received. 


(c) X-Ray Methods 


For the alloys with titanium, vanadium, chromium and cobalt, 
filings were prepared and strain annealed at 600°-800°c, whilst for 
the iron—manganese and iron—nickel alloys the treatments were at 
460°—700°c in order to avoid the formation of 2-phase (x+y) structures. 
The filings of alloys with titanium, vanadium and manganese reacted 
with silica and were therefore placed in alumina protectors inside silica 
tubes for the strain annealing treatment, and precautions taken to 
remove adsorbed gas before sealing off. Owing to the ferromagnetic 
nature of the alloys possible fragments of the file could not be removed 
magnetically, but as the alloys were soft it is improbable that appreciable 
errors were introduced in this way, and the work of Owen and Liu (1949) 
has shown that filings and solid specimens of iron-nickel alloys give 
identical lattice spacings. 

The X-ray specimens were prepared in the form of thin cylinders, using 
Canada Balsam as adhesive. The powder photographs were taken in two 
19cm Unicam cameras which were checked against one another, the 


‘a 
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exact lattice spacings being determined by standard extrapolation 
methods. The coefficient of expansion of iron was used to correct the 
values to refer to 20°c and all results are given in kx units, uncorrected 
for refraction, and subject to this condition should be accurate to better 
than +0-0001 kx so far as x-ray technique is concerned. The standard 
wavelengths on which the present results are based are : 


Cobalt Ka, A=1-78529 kx, 
Ko. A=1-78919 kx. 


Table 1. Results for Iron Alloys at 20°c 


: Lattice 
Solute Atomic ° 
Bone solute a tle eter Notes 
x 


Pure Fe 2-8604(3) Mean of seven determinations 


Arc-melted alloys 


VAro-melted alloys 


N.P.L. alloys 


Induction-melted alloy 
N.P.L. alloy 


\nduction-melted alloys 


N.P.L. alloy 
Induction-melted alloy 
N.P.L. alloy 


SHwWodTSSH 
SARNOS HE 


Induction-melted alloys 


OS & bo COROT ee 


ae Oe 
Dwr © 


Master alloy 


N.P.L. alloys 
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$3. Tue Larrice Spactne or Pure [Ron 


The lattice spacing of the pure H iron was determined as 2-8604(3) kx 
at 20°C in good agreement with values of between 2-8601 kx and 
2-8605 kx obtained by different previous workers after reduction to 20°c. 


$4. Larrick SpacInGs oF Sotrp SoLuTIONS IN [RON 


The lattice spacings obtained in the present work are summarized in 
table 1 and the following points may be noted. 


(a) Iron—Titanium 
No previous accurate data existed for this system, and the present 


results are shown in fig. 1. Of the eight points, seven lie on a straight 
line, and one is incorrect for reasons which are not understood. The 


Fig. 1 


Fe-Ti 


At. % Ti 


straight line does not, however, extrapolate to the value for pure iron, 
and the first additions of titanium produce a relatively small lattice 
expansion. ‘This effect is not understood, and is too great to be due to 
intercrystalline adsorption of solute. It is also too great to be ascribed 
to oxidation or nitriding during melting which might convert titanium 
into oxide or nitride remaining in a form which was dissolved during 
chemical analysis and so give an apparent titanium content higher than 
that of the solid solution. The effect resembles that found by Axon 
and Hume-Rothery (1948) for solid solutions of high purity (99-97%) 
magnesium in aluminium, and here later work by Poole and Axon (1951) 
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showed that the effect was removed by using more pure (99-98°) 
magnesium. It is probable, therefore, that the effect is due to traces of 
impurity in the alloys, but in spite of much work it is not understood. 

In fig. 1 the dotted line shows the results expected if Vegard’s law were 
obeyed, and the lattice spacing of titanium were that of the P (body- 
centred cubic) modification corrected to refer to room temperature. In 
the iron alloys there is a contraction compared with Vegard’s relation, 
and there would still be a slight contraction if the straight line part of 
the present curve were displaced to pass through the value for pure iron. 


(6) Iron—Vanadium 
In fig. 2, the full points show the lattice spacings of the two arc-melted 
alloys, and these lie on a straight line passing through the value for pure 
iron. The open points refer to alloys which were accidentally very 


slightly contaminated by tantalum and which are therefore not included 
in table 1. The dotted line shows the lattice spacings expected according 
to Vegard’s law, and there is again a marked contraction. 


(c) Iron—Chromium 
Figure 3 (a) shows the results obtained for iron-chromium alloys, 
whilst the dotted line shows the lattice spacings expected according to 
Vegard’s Law. Figure 3 (b) shows the present results compared with 
those found by Preston (1932). When allowance is made for the slightly 
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different values for pure iron, there is good agreement between the two 
determinations although the present results indicate a slightly smaller 
increase in lattice spacing. In this system the lattice spacings are greater 
than would be expected from Vegard’s law. 


Fig. 3 

2.864 reas 2-864 
2.863 kX 
kX 
2-862 2.862 
2-861 

——®— present authors 

— & — Preston 
2.860 2.860 

At.% chromium 
(a) (b) 


(d) Lron—Manganese 

Figure 4 (a) shows the results obtained for iron-manganese alloys, 
together with those of Troiano and McGuire (1943) (fig. 4(b)). The 
points for the seven alloys examined in the present work lie on a straight 
line to a very high degree of accuracy, but this straight line extrapolates 
to a point about 0-0004 kx below the value for pure iron, although the 
effect is not so pronounced as for the iron-titanium alloys. If allowance 
is made for the slight difference in the lattice spacings of the iron used, 
the present results agree well with those of Troiano and McGuire which 
also show an abnormality at low concentrations. The data suggest a 
slight contraction compared with those expected from Vegard’s law, shown 
by the dotted line in fig. 4(a) although this difference would almost 
disappear if the straight line portion of the curve were displaced parallel 
to itself so as to pass through the value for pure iron. 


(e) Lron—Cobalt 
Figure 5 (a) shows the results obtained for iron-cobalt alloys together 
with those of Rodgers and Maddocks (1939), and of Ellis and Greiner 
(1941) (fig. 5(b)). The lattice distortion is very small, but is greater 
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than that found by Ellis and Greiner, although in good agreement with 
that of Rodgers and Maddocks if allowance is made for the difference 
between the lattice spacings of the iron used. In fig. 5 (a) the upper 


Fig. 4 


Fe-Mn 2.862 


kX 
2.86) 
—e— present authors 
— #- - Troiano a McGuire 
2-860 
O { 2 3 4 
At.% manganese 
(a) (0) 
Fig. 5 
Fe-Co 
2.8610 
a) 
2.8608 a _A o 
Ee ie are 
A ; 
2.8606 ee 
B60} kX A 
2.8604 z 
Nout 
oe 
2.8602 a 
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dotted line shows the results expected if Vegard’s law held, and the 
atomic diameter of cobalt were that observed in «- and f-cobalt. This, 
however, refers to structures with coordination no. 12, and should 
therefore be reduced by 2-5°/,* in order to give the values appropriate 
to a hypothetical body-centred cubic modification of cobalt. If this 
correction is made, the results show an expansion compared with Vegard’s 
relation (lower dotted line, fig. 5 (a)), although there is a slight contraction 
if no correction is made. 
(f) Lron—Nickel 

Figure 6 (a) shows the results obtained for iron—nickel alloys, together 

with those of Bradley, Jay and Taylor (1937), and of Owen and Sully 


(1941) (fig. 6(6)). When allowance is made for the different values for 
pure iron, the present results show a smaller expansion of the lattice 


Fig. 6 


Fe-Ni 


2.864 


kX 
2.862 

—e— present authors 

- #&- Bradley, Jay « Taylor 
2.860 


O° \ 2 3 4 5 
At.% nickel 


(a) (b) 


than was claimed by Bradley and Jay, but are in good agreement with 
the data of Owen and Sully. The lattice spacing/composition curve is 
at first a straight line, and then shows a decrease in the rate of expansion. 
The initial straight portion of the curve shows an appreciably greater 
rate of expansion than would be expected from Vegard’s law for a 
hypothetical body-centred cubic nickel with an interatomic distance 
equal to that of face-centred cubic nickel (upper dotted line, fig. 6 (a)), 


whilst if the interatomic distance in nickel is corrected to allow for the 
eis nT eRe ie ie Te a 
* See p. 1304 for details. 
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contraction on change of coordination no. 12->8, the positive deviation 
from Vegard’s law is even greater (lower dotted line, fig. 6 (a)). 


§ 5. Discussion 
Figure 7 shows the lattice spacing/composition curves obtained in the 


present work, together with that of Norton (A.S.M. Handbook) for the 
restricted solid solution of copper in «-iron. Considered empirically, the 


steer) 


2.867 Ti 


2.866 


kX Cr 


2-862 Mn 


2.861 .. 
Lf Co 


2 3 4 5 6 a 
At.% alloying element in a-iron 


curves fall into two sequences, in which at a given atomic percentage of 
solute, the lattice expansions are in the order Ti> V>Cr and Co< Ni<Cu. 
The behaviour of Mn is anomalous, its effect being nearly the same as 
that of chromium.” 

Figure 8 shows the lattice spacing/composition curves of the same 
elements dissolved in nickel, and here the effects are in the same general 
order Ti>V>Cr and Co<Cu, although the lattice of nickel is very 
slightly expanded by cobalt, in contrast to the order Co<Ni<Cu found 
in the iron alloys. For the nickel alloys, the abnormal behaviour of 
manganese is more pronounced, and the Ni-Mn curve lies close to that 
for Ni-Ti. 
ano Le ES eS eS es See ee ee eee) rae 

* Detailed discussion of the curve for manganese is difficult because of the 
abnormality at low concentrations. If the Fe-Mn curve were shifted parallel 
to itself so as to pass through the value for pure iron it would lie between the 
curves for Fe-V and Fe-Cr. 
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It is also instructive to compare the data for the iron alloys with those 
found by Hellawell and Hume-Rothery (1954) for solid solutions of the 
corresponding elements of the Second Long Period in the close-packed 
hexagonal solvent ruthenium. Here at a given atomic composition, 


Fig. 8 


3.56, 
Mn 


3.55 ip 


Gui 
Fe 


Co 


3:5) 
12) 5 [@) IS 


At. % alloying element in nickel 


the expansions of the ¢ lattice spacing of ruthenium are in the order 
Zr >Nb>Mo, and Rh<Pd, the effect of Rh being small,* and that of 
Pd slightly less than that of Mo. Quite different effects are shown by 
the a spacings, which are contracted by Zr, and very slightly contracted 
by Rh, whilst Nb and Mo produce large and roughly equal expansions, 
and Pd a less marked but considerable expansion. Considered empirically, 
therefore, it is the effects for the c spacing of ruthenium which resemble 
those found in the present work for solid solutions in iron.+ 

To discuss the curves for iron and nickel in greater detail, correction 
must be made for the effects of coordination numbers. The data for 
iron (Basinski, Hume-Rothery and Sutton 1955), and manganese 
(Basinski and Christian 1954) suggest that Goldschmidt’s empirical 


* Data for Te alloys are lacking. 

+ In the ruthenium alloys the effects of the solutes upon the volume of the 
unit cell are such that Zr causes a contraction, Rh has no effect, whilst the 
remaining elements cause expansions in the order Nb, Mo (approx. equal) > Pd. 
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factor of 3°% contraction in atomic diameter on passing from face- 
centred cubic to body-centred cubic coordination is slightly too great, 
and that a contraction of 2:5°% is in better agreement with the facts. 
Using this correction, the values in table 2 are obtained for the atomic 
diameters of the elements concerned. 


Table 2 


Klement 4h V Cr Mn | Fe Co Ni Cu 


Atomic diameter in | 2°89" | 2-704 2-56f | 2-61 | 2-54] 2-50 | 2-49 | 2-55 


face-centred cubic 9.94 


Atomic diameter in | 2:84} | 2-62 | 2-49 | 2-54 | 2-48 | 2-44} | 2-43 | 2-491 
body-centred cubic 


-* Values for the two close distances of approach in the close-packed hexagonal 
structure of a-Ti. 
+ Value from extrapolation of lattice spacings of solid solutions in B-Ti. 
+ Values calculated by assuming 2-5%, contraction on passing from face- 
centred cubic to body-centred cubic. 


Using these values, the data for Fe-Ti and Fe—-V show a contraction. 
compared with Vegard’s law, whilst those for Fe—Cr, Fe-Co, Fe—Ni and 
Fe—Cu show an expansion. This difference is reflected in the equilibrium 
diagrams of the systems concerned. In the systems Fe-Ti and Fe—-V 
there are intermetallic compounds or intermediate phases of high melting 
point, which suggest strong bonding of the two kinds of atom. In the 
system Fe-Cr, the equiatomic o phase is stable only below 815°c, whilst 
in the remaining systems there are no intermediate phases apart from 
superlattices stable only at low temperatures. 

Although there are departures from Vegard’s law, the lattice expansions 
Ti>V~>Cr are in the same order as that expected from the atomic 
diameters, but from this stage onwards the order is quite different, 
and the behaviour of manganese is clearly anomalous. For the elements 
succeeding iron, the lattice expansions are in a well defined sequence 
Co<Ni<Cu, and bear no relation to the atomic diameters which are 
all nearly the same, and to a higher degree of accuracy are in the order 


Ni 2:43<Co 2:44<Fe 2:48<Cu 2:49. 


We may use the term “outer electrons’ to denote all the electrons 
outside the argon shell, and ‘ electron number ’ to describe the number 
of electrons per atom concerned. in metallic bonding.§ In the ‘ collective 


§ This has usually been called the valency of the atom in the metal concerned, 
but the term has been used with so many meanings that we prefer to speak of 


the electron number. 
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electron theories’ many authors now regard all of the outer electrons as 
taking part in the bonding process, so that the electron number is equal 
to the number of outer electrons. From this viewpoint, the break in 
the sequence of physical properties in the middle of the transition elements 
is the result of an increase in the repulsive part of the wave functions 
of the d-electrons, which becomes very marked on passing through 
Group VIII (Co+Ni, Rh-+Pd, Ir+Pt). In the resonating bond theories, 
the cohesion in crystals of the transition metals results from the resonance 
of hybrid spd bonds, and the breaks in the physical property sequences 
are the result of some electrons entering atomic or non-bonding orbitals. 
In the Pauling (1938) hypothesis the electron number is regarded as 
increasing from 1 to 5 on passing from KV, and then as remaining 
constant at 5-78 on passing from Cr+Ni. As poimted out by Hume- 
Rothery, Irving and Williams (1951) this scheme is not in agreement 
with the physical property sequences, and these authors suggest that, 
in the First Long Period,* the electron numbers increase from 1 to 4 
on passing from K-+Ti, with an uncertain value for vanadium in Group VA, 
and a value for chromium which is lower than those for molybdenum and 
tungsten in the later Periods. The abnormal properties of manganese 
were regarded as associated with more than one electronic state, whilst 
the relatively large interatomic distance in y-Mn and the low melting 
point of 6-Mn were associated with an electron number 2, corresponding 
with the very stable divalent Mn?+ compounds. Later work by Basinski 
and Christian (1954) has shown that 6-Mn has also a relatively large 
closest interatomic distance. All the physical properties suggest a rise 
in electron number on passing from Mn-Fe, together with a very gradual 
decrease in the sequence Fe+Co-+Ni. Hume-Rothery, Irving and 
Williams pointed out that the electron numbers suggested by the physical 
properties ran parallel to the predominant valencies in the chemistry of 
these elements. 

In an early discussion Mott (1944) regarded copper, silver and gold 
as consisting of ‘ hard ions ’ pulled into contact by the s valency electrons, 
with a simultaneous additional cohesion of the van der Waals type 
arising from the perturbation of the outer d electrons of the ion, and 
this view has been discussed by one of us (Hume-Rothery 1955 a). 
Recently Mott (1953) has suggested that this concept can be applied to 
metals at the end of the transition series, the van der Waals perturbation 
being treated as an excitation of the d—p-type. This picture, although 
not its mathematical details, is not so very different from that of the 
resonating bond theories, since each imagines a core of repulsive electrons, 
and each involves s, p and d states for the cohesion. On passing back 
from Ni+Co-+Fe the d electron clouds become larger and more loosely 


* In the Second and Third Long Periods the electron numbers increase from 
| to 6 on passing from Groups LA->VIA, remain high in Groups VII and VILLA, 
and then fall rapidly on passing to Groups VUIIB and VILIC. 
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held, and the van der Waals perturbation forces increase until they 
become so large that they break into normal metallic cohesion. 

When we try to combine the different concepts, we may regard the 
almost exactly equal atomic diameters in the elements Fe, Co, Ni, Cu 
as the result of the following opposing tendencies : 

(1) With increasing atomic number the 3d electrons become more 
firmly held, and the electron clouds of both 3d and 4s electrons contract, 
and thus tend to produce a smaller atomic diameter. 

(2) The contraction and firmer binding of the 3d electrons means that 
they are less easily perturbed, and this reduces the attractive force, 
whether the latter is regarded as resulting from spd hybridization, or 
from van der Waals perturbation of the d electrons. At the same time 
the filling up of the (3d) shell increases the exchange repulsion for a 
given degree of overlap of the electron clouds, and these two effects 
tend to increase the atomic diameter. 

(3) The fall in electron number suggested by the chemical valencies 
on passing from Fe-+Cu weakens the attraction. In so far as this fall 
results from the withdrawal of attractive d electrons it is perhaps really 
included in (1) and (2) above. 

Using the low electron numbers of H-R, I and W, a considerable 
proportion of the d electrons round each atom in an iron crystal is 
repulsive or non-bonding, but the number of these repulsive electrons 
(5 to 6 per atom) is only slightly more than one-half of that required 
to fill the (3d)1° group, and so considerable overlap of the cores can take 
place. The atomic diameter of copper when corrected for coordination 
number is almost equal to that of iron (table 2) so that if we regard 
copper as consisting of hard ions with completed (3d)!° groups in contact, 
the insertion of a copper atom into the iron lattice will produce an overlap 
of the completely filled (3d)1° electron cloud with the (3d) electron cloud 
of the iron atoms, and the resulting exchange repulsion may be regarded 
as the cause of the marked expansion of the iron lattice, this effect being 
reinforced by a weakening of the van der Waals attraction,* and a 
decrease in the electron number of the solute (Cu) as compared with the 
solvent atom (Fe). With nickel as solute we may expect the same 
effects but to a lesser degree. The (3d) core of a nickel atom in solution 
in iron is probably nearly although not completely} filled, and although 
the unfilled states permit bonding, any marked overlapping will result 
in exchange repulsion. With cobalt as solute the (3d) core is less 
completely filled, and exchange repulsion will be less marked. The 
curves in fig. 7 show clearly that although the expansions produced by 
a given atomic percentage of the solutes are in the order Co<Ni<Cu, 


* Or alternatively a decrease in the perturbation which enables metallic 


bonding to occur. 
+ It is possible that the holes in the d shell of the nickel atoms are filled by 


electrons from the iron atoms (see p. 1308), 
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the effect of cobalt is very small, and we may conclude that the exchange 
repulsion increases enormously when nickel is reached. The fact that 
cobalt produces a slight lattice expansion of iron, corresponding to a 
marked positive deviation from Vegard’s law, can readily be understood 
as partly due to the exchange repulsion effect, and partly due to the 
factors outlined in (2) above. The atomic diameter of iron (2:48 kx) 
results from an electron number slightly greater than that of cobalt, 
reinforced by the van der Waa!s attraction of the iron d shell core, 
which is larger and more easily perturbed than the d shell core of cobalt. 
The substitution of cobalt for iron thus weakens the van der Waals 
attraction and increases the exchange repulsion. 

If the above interpretation is correct, we can understand why the 
lattice expansions of iron produced by Co, Ni and Cu bear no relation 
to the atomic diameters of the solutes. We can also perhaps understand 
why the lattice of nickel is slightly expanded by cobalt, and more so 
by iron, because in nickel the (3d) shell is so nearly filled that exchange 
repulsion may predominate over increased attraction when the larger (3d) 
electron clouds of cobalt or iron are substituted.* The same interpretation 
enables us to understand why the lattice of nickel is expanded by copper 
more than would be expected from Vegard’s law, because, although the 
3d electron cloud of a copper atom is smaller than that of a nickel atom, 
the fact that the (3d)!° sub-group is filled in copper means that intense 
repulsion will result from a relatively small degree of overlap, particularly 
since the magnetic evidence shows that the (4s) electrons of the copper 
atoms fill the holes in the (3d) shells of the nickel atoms and hence make 
the latter more liable to give rise to exchange repulsion. When nickel 
is dissolved in copper, there is a positive deviation from Vegard’s law 
which can be understood as the result of exchange repulsion, but the 
lattice is contracted and not expanded. It is perhaps significant that, 
as emphasized by Coles (1952), nickel when dissolved in copper gives 
paramagnetic alloys, indicating that the (3d) shells of the nickel atoms 
are not being filled by electrons from the copper atoms. 

From the above viewpoint, the distortion of the iron lattice produced 
by manganese is anomalous, and the effects of cobalt, nickel and 
copper result from the interplay of several factors. The effects of 
titanium, vanadium and chromium are more simple, and of these three 
elements the size factor is really favourable only for vanadium and 
chromium. In recent work in this laboratory by A. Hellawell, the 
liquidus and solidus curves of this series of iron-rich alloys have been 


ee a 

* The magnetic evidence shows that the saturation moment of nickel is 
increased by solution of cobalt or iron, but it is not possible to say whether 
this is accompanied by a filling of the holes in the 3d band of the nickel atoms. 
In view of the strong tendency to fill the d shell in the last element of each 
transition series, it is possible that the substitution of an earlier element 
(Co, Fe, Mn...) always tends to fill the 3d shell of nickel, and so to increase 
the exchange repulsion. In this case the increase in saturation moment would 
result from the solute atoms, 
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determined accurately, and the results suggest that in dilute solutions 
a given atomic percentage of vanadium depresses the melting point by 
slightly less than twice* the depression caused by the same atomic 
percentage of chromium. The present work shows that the lattice 
expansion for a given atomic percentage of vanadium is slightly less 
than twice that for the same atomic percentage of chromium. For these 
two systems, therefore, the lattice expansions and depressions of melting 
point are related by the same factors, as is the case for some silver alloys 
examined by one of us (Hume-Rothery 1955 b). This correspondence, 
however, no longer holds for the system iron—titanium where the lattice 
distortion is much greater. 
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SUMMARY 


The eigenvalue problem which describes the onset of convection in a 
fluid sphere heated within is solved by means of expansions of the 
convective velocity field and temperature in spherical Bessel functions. 
The Rayleigh number C,,,,, at which there first appears the nth radial 
convection pattern with the spherical symmetry of the /th spherical 
harmonic, can be computed for any / and n by this device, and bracketed 
to-arbitrary accuracy. A table of C,,,, for / from 1 through 6 and n=1 
and 2 is compared with Chandrasekhar’s values of C 1,1 obtained by a 
variational method. A similar attack is shown to succeed on the classical 
Rayleigh—Benard problem of convection in a horizontal layer of fluid. 
Finally, the class of problems which aids be amenable to such an 
attack is delineated. 


§ 1. INTRODUCTION 


THE linear partial differential equations which describe the onset of 
convection in a fluid sphere heated within were first derived by 
Wasiutinsky (1946). These equations constitute an eigenvalue problem, 
the eigenvalues being the Rayleigh numbers at which various patterns 
of convection can set in. Jeffreys and Bland (1951) computed the 
lowest eigenvalues by a variational procedure, and by means of a more 
incisive variational technique Chandrasekhar (1952; this paper will be 
referred to hereafter as I) was able to obtain quite accurate estimates 
of the lower eigenvalues. Variational estimates of eigenvalues suffer 
from two disadvantages: they are upper bounds, which in principle 
cannot bracket the true values; and the higher eigenvalues are 
inaccessible without inordinate accuracy in the computation of the 
lower ones. Therefore, in this paper the eigenvalues will be computed 
by a more laborious technique which does bracket them and which can 
be used to compute any eigenvalue independently of the others. This 
method consists of an application of Chandrasekhar’s scheme for lowering 
the order of a high order differential equation. The method exhibits the 
eigenvalues as the zeros of a certain explicit meromorphic function of a 
single complex variable and, as a by-product, yields expansions for the 


* Communicated by S, Chandrasekhar, F.R.S, 
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eigenfunctions (the convection currents) in series of spherical Bessel 
functions. Chandrasekhar’s formulation of the problem will be used 
throughout this paper, and with two exceptions his notation will be 
followed in detail. Just as his procedure can be extended to spherical 
shells (Chandrasekhar 1953 b), so can the techniques to be presented 
here. However, this extension involves laborious computations, so it 
will be omitted. 


§ 2. THE EQUATIONS GOVERNING MARGINAL STABILITY 


Chandrasekhar’s derivation of the linearized equations for the onset of 
convection in a fluid sphere heated from within appears in I and will 
not be reproduced here. In copying his equations, their numbers in I 
will be given beside their numbers in this paper. In his notation, x is 
the thermometric conductivity of the fluid, v is its kinematic viscosity, 
p is its average density, « is its volume coefficient of thermal expansion, 
and 6x8 is the heat generated per unit volume per unit time in the fluid. 
The universal constant of gravitation is G, and y=(4/3)7paG, while 
L=—i(rxv) is the ‘angular momentum’ differential operator which 
arises in connection with the spherical harmonics. Finally, uw, is the 
radial velocity of the fluid, and @ is the deviation of the fluid’s 
temperature from the steady temperature in the absence of convection. 
Chandrasekhar’s equations, based on the principle of exchange of 
stabilities, are then 


eV2G=—26ru,, . . . « (1) (1 (24)) 

and 
V4(ru,,)=(y/v)L?0. ee (2) ae (34) 
Taking R to be the radius of the sphere, the boundary conditions are 
G=20 av r= f, ee ee Ore (>). ((16)) 
fi =) ar, ete ean an mee (4)) 6 (17(1.6)) 
and. £ (ru,)=0 at r=R for a rigid boundary, . (5) (I (19)) 
or a (ru,,)=0 at r=R for a free boundary. . (6) (I (23)) 


or 
These equations can be separated by the substitutions 
ru,= Wr) Y,"(3,¢) and 0=O(r) Y"(3, ¢). 


The resulting differential equations are 


1p@ +1) 2BW 
? Ted ia 2 ee 
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When /=0, the condition 

vV.u=0 - a Se 
becomes d(Wr)/dr=0, which implies that Wr=r?u, must be a constant. 
In the absence of sources at the origin this constant must be 0, so there 
will be no eigenvalues unless />1. 


§ 3. THe SOLUTION OF THE EIGENVALUE PROBLEM 
Take *=r/R as a new independent variable and immediately omit the 
tilde. Define 


wg l(d+-1) 
cae ory te «: A) 439 » <n 
2 
and Ch BY Ro. 5 a 4 aoe! i 
KV 
Then (7) and (8) become 
1p d? Ul+1) : 
Aon | rF=—Ul+1)C, »W, ). 29h ae 
Lfedtas Tle l)iaee 
and =| — al A ed os +e OSes 
The boundary conditions, obtained from (4), (5) and (6), are 
Wit)h=F(1)=0,- . . .. . . « (14) ee 
W and F behave like rv’; asr+0, . . (15) (I (44)) 
and W’'(1)=0 for a rigid boundary, 
be, (16) (I (47)) 
or W’’(1)=0 for a free boundary. 


Equations (12) and (13) together with boundary conditions (14), (15) 
and (16) constitute the mathematical problem to be discussed. The 
eigenvalue C,,,, is called C, in I because Chandrasekhar is interested 
only in the lowest characteristic value for each /, denoted here by C, ,. 
Since the spherical Bessel functions j,(z)=(7/2z)"?J,,, 4/9(z) satisfy 


Ly ah Ker Ln 
ae se | silo) = —a2jer) ols an 


they form a convenient basis in terms of which to expand F and W. In 
fact, as Chandrasekhar has pointed out privately to the author, if ® 
and ¥ denote real and imaginary parts of a complex number, then 
W=Aj,(ar)+ BRL j,(ar e)+(—1)ij(ar e) 
+ OBL jy(ar e9)—(— 1h, (ar o*)] 
where A, B, C are real numbers chosen to satisfy (14), (15) and (16), 
and « is any real root of the characteristic equation determined by (14), 
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(15) and (16). For example, if the boundary is rigid, this equation is 


Jul) Bj ,(a e'/3) + (— Lyju(ot ei27/3)) 
Dl jla eo?) —(—1)49 (a e*)] 
ga | Jie) «=»: Ble (x 0°?) -(— 1) e789 (a 0") 
Hlett*/3 (a e?7/3) _ (— 1) ef27/35 (a e'27/3)) 
Jriarla) Ble), (a e*/9)+ (— 1) e797, (x 07/9)] 
Her}, (a. e#7/3)—(—1) e734, | (a e/8)] 


(18) 


Then C,,,—«°. Equation (18) is sufficiently formidable to prompt other 
approaches. 


The Sturm—Liouville theory ensures that / can be expanded as 
tee) (Oy, ae eee ae ae, ace a. (LD) 
n=1 


where «, ,, is the nth positive root of j,(2)=0. We shall usually omit the J, 
writing simply «,. The boundary condition that F behaves like r’ as 
r—+0O ensures that (19) converges uniformly on the interval [0, 1j 
(Watson 1944, p.615). Because of this uniform convergence, the first 
to fourth anti-derivatives of F can be computed by integrating (19) term 
by term. Therefore eqn. (11) together with the boundary conditions (14) 
and (15) can be solved exactly for W in terms of F’, by virtue of (17). 
Thus the sixth order problem embodied in (10) and (11), namely 


1p d? JWi+1)78 
=| 7a S | rW=—UI+1)C,,W, . .-- ~ (20) 


can be reduced by four exact integrations to a second order problem. 
This technique was introduced by Chandrasekhar, who discusses many 
of its applications (see particularly Chandrasekhar 1954). In I he applies 
the method to eqn. (20) as outlined above. It should be observed that 
either F or W can be expanded in a uniformly convergent series of spherical 
Bessel functions as in (19). We expand F because then the boundary 
condition remaining after the four exact integrations is the one that takes 
two possible forms, namely (16), and the necessity for discussing two cases 
is delayed as long as possible. Also, starting with F facilitates comparison 
with Chandrasekhar’s procedure in I. 

The result of solving (11) for W in terms of the F given in (19) is 
Wey EAM jfoqr) BA EOME Te ta. 21) CL ) 
Note that the 4, in (16) differs from Chandrasekhar’s 4, in I by a factor 
(2/m«,,)/2. Boundary condition (13) requires M=N=0, and W(1)=0 
implies B= —C, so that (21) becomes 


Woj= EF Aaj fayr)+ Br)... (22) (173) 


n=l Sn 
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The requirement F(1)=0 now implies 


q < Ze © . 
Be 4 a 5 ad i+ 1%n)> . 2. (23) (iia 
n=1 “n 
where qg=2 for a rigid boundary, 
and q=—4/(21+1) for a free boundary. 


Now (22) can be substituted in (12) and, since (22) converges uniformly 
if (17) does so, the result can be integrated term by term to produce 


ae 2) A me yl t4 ylt+2 : E 
Ne 2, eho |garsa— garry +P 
The eae condition (15) makes H=0, and (14) determines D, so that 

Et) gape ae pita ylte (20--7)r? 
en, 2 z,soient) TB Fees ~ 2(21+3) * 4(21+-3)(21-+45) 
(24) 


Now (19) and (15) imply 


2 1 i} : 
a Aa gerTe me =) Jin?) 


n= i 


= yl +4 plt2 (20-+-7)r? 
= Fes; ~ 3@r3) 7 Tas | - (25) 


Since the Fourier—Bessel expansion of the right side of (25) is easily 
obtained, (25) determines 4, in terms of «,, B and C,,,,. 

To obtain A, in terms of the Fourier—Bessel coefficients of the right 
side of (25), we use the following standard relations among spherical 
Bessel functions : 


el 


J elena ml 2dr= sia (07) Ovens D ° ° s ‘ , © < (26) 
ey 

| miler) dr=(on\jrsr(n)s © 2 2 2 2 ee we es QM) 
ri ; ) Se 2(27+-3 

[Slay dr— ean] AEE =|. Oy en re 


and 


2 (29) 


1 nn n 


i: Hj aq) dr — Peay ital eee) 7 
my 0 v 


Multiply (25) by 7?j,(«,7) and integrate from 0 to 1. Using (26) through 
(29), the result is 

4a,,(21+-3)l (I+1)C l,m 
Jrgr(%y Lx n — (i+1)C tbs mle 
If C,,, is known, the corresponding temperature © and radial velocity W/r 
can be computed from (30), (19), (22) and (23). All the results of this 


od eee (30) 
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paper follow from (30). In particular if (30) is substituted in (23) and 
B is cancelled, the following equation for C, ,, is obtained : 


“lm 


l - i-+-1) De 


TeTaca\ —— mews (81) 
Q(21+-3) pay Mn LUE IVC), m— On 
§ 4. GENERAL REMARKS ABOUT THE LOCATION OF THE HIGENVALUES 


The remainder of this paper, except for §§ 7 and 8, is a discussion of 
the roots of the equation 


i 
where b(z)= 5 cde dena SSRN ITD dla (33) 


n= a) also sy) 


Figure 1 is a graph of #,(z) for real z. If 7,,,, is the mth positive root of 
(32), then 

Hel Ge a =r peer wt", oat (34) 
The subscript / will be omitted from 7,,,,, hereafter. 

If the interiors of small circles around the simple poles z=z,,° are 
deleted, then in any bounded region of the complex plane the series (33) 
converges uniformly. Therefore (Copson 1935, p. 97) its sum is analytic 
and can be differentiated term by term. Thence 

o YL. 4 


n=1 (230, 


It is clear that if z is real /’(z)<0, so (31) has precisely one real root 
between « y® and «y,,%. Furthermore, if <0, %(z)>0 and #(z) decreases 
monotonically from 


Ray ‘ 
1m ye aT 
2 6 
2>—-0 n=1 %n (20, ) 
: z 
Since lim 5 =1 
Z—H 
Z>-—- oO ~ n 


exists uniformly in n, the former limit can be computed term by term 
(Knopp 1928, p. 338). Therefore 


ae L 
lim d&(z)= 2 sei 
Z>— 0 el Hi 
But (Watson 1944, p. 502) 
oD 1 1 
= FT oT Sa) ° 36 
or ey 2(21-+-3) ( ) 
il 
lim 1 reer a SAA ey ae Se ee) 
eae jim h@)= serp3) 


The value of this limit is precisely the value of g(2/+1) for the rigid 
boundary problem. Therefore (35) and (37) together imply that in both 
the free and the rigid boundary problems, all real roots of (31) are positive. 


1316 G. E. Backus on the 


Since 
: ©  r(r—a,,° COS 8) ae - es 
wre”*)= LY ———“*.—_, at Tea sn as , (38) 
n=1 %% |r x, © | n=1 | T— Oy @ 


eqn. (32) can have only real roots. It is thus possible to conclude that 
all the roots of (32) and hence all the eigenvalues (,,,, are positive real 
numbers. In the free boundary problem there is one extra root of (32), 
between 0 and «,°, to which there corresponds no root of the rigid 
boundary problem. In what follows, the root of the free boundary 
equation between a) y and a) y,, will be called 7, y,,. The root of the 
rigid boundary equation which lies between xy and ay, will be 
called r, y. Thus in both problems the roots are 7, 1<7, 9<?7j,3<.... 


Fig. 1 


1 
Freyery) FOR] RIGID BOUNDARY 


Or) Sites | ee ee 
0-0 
-0.1 
Gata a) FOR| FREE BOUNDARY 
-0:2 
-0.3 
=2 Sh @ |) 2 4 6 8 10 12 14 6 7's 


It is clear from what has been said that for a fixed /, if7, ,, is approximated 
by x, ,,°, the percentage error thus committed approaches 0 as n approaches 
co. Watson (1944, p. 506) shows that «,,, is asymptotically equal to 
(n-+-l/2+-,)m for large n, where A, is an integer which depends only on J. 
Thus for large n, 7, ,, is asymptotically equal to (n+-1/2+-d,)®r. 

Another obvious question is the behaviour of r,,, for fixed n as +o. 
We shall consider only 7, ,. Chandrasekhar shows that 


= ay VI 41/9(% 1) : , 
11S OT eo 
Be “a 412%, DQ 45 /2(%, 1) ( ) ( ( )) 


By means of the well-known recursion relations among Bessel functions 
(Watson 1944, p. 45) this inequality can be rewritten 


(1-S)<-, Sie) eee 


6 2 re 
O14 4 red 
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For the rigid boundary problem 7, ,>«, ,°, so 


5) ‘ 
! 0-2 )\<-< ie Pe OL) 
} 


6 2 z = 6 
Hy Hy ait CSA 


for rigid boundaries. 
Watson (1944, p. 485) shows that 


Dee teeny veo ce isl Sin whian, (22) 
Then (41) implies 
2(21+3) a, 1° 
1— ——, < >< l Po EM eee eS 
(t-r4)? - Liat Pe ce 
for rigid boundaries. This means that as /> 0, 1; 1s asymptotically equal 
to x, ,° and therefore to (J+ 3)° (Watson 1944, p. 521). 
For the free boundary problem more care is required. From (33), 


Deer hese 280 (44) 
a n=1 0,,°(1—z/a.,,°) 
If [2 | <a, 
ce) 2 oO ~ k nD 0 ak 
i Dal 5 Oho BE AS 
a ee Cn et (=) a hoa oo tre ( 5) 


The Weierstrass double sum theorem (Copson 1935, p. 97) then ensures 
that the order of summation can be reversed, so, if |z|<21°, 


UO a ara a ee eee (46) 
k= 
where 
iva) | 
o, : (47) 
k ae era 
Clearly, 


, es NN O42 
—p(z)<oy 2 (=) = = =. alto (48) 


K=0 \%1 
Therefore the lowest root, 7, ,, of (32) is larger than the root z of 


Cs (21+ 1) (49) 
ne © MER) Nagler 

ree 

A(21By oe hehe £50) 
(14 (I) Ory n) 


The root of (49) is pe 


Combining this with (40), 
1 4(21+-3) 1 1 ( (ioe 
ale Raho nee ae eee ie (See 51) 
ny (14 (20-1 )o;, fs * Ty t m1 i° wi 4(21-+-3) yee ( 
for the free boundary problem. Now (51) implies that for the free boundary 
problem 1, ; 1s asymptotically equal to %, 1° as > co, with a per cent 
error which behaves like 1/1 as l>o. To see this, use the fact that 
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is a decreasing function of n (Watson 1944, p. 515) which 


on+1 &n 
approaches 7 as ” approaches oo (Watson 1944, p. 506). Then 
cal gee | 1 x 1 = 1 iS fs dx 
ae vitae cece mee RS een SSS SS, 
‘es ey re ay! ¢ n= (%1+n7)* ~ 48 w/z (%1-+7@)° 
1 1 
~ ae Fale) 
Hence Bae ay” meld 1 fee 
Tt Si | vag ae 2 (B21 


Since «, ,>/+-4, the assertion is established. 

One remaining general remark about the distribution of the eigenvalues 
is that among all eigenvalues of the problem specified by eqns. (1) through 
(5) or (6), the smallest is C, ,. This assertion is of some interest since the 
lowest of the C,,, is presumably the Rayleigh number at which convection 
will begin in an initially cool sphere which begins to warm up internally ; 
and the lowest C,,,, is also the Rayleigh number at which convection will 
cease in a sphere of fluid losing heat through its surface. 

Clearly Cy m<Cy,m41 SINCE Ty, m<l1,m412 and by eqn. (34) Cy, m=",m/Ul+ 1). 
There remains only to show that C,,<C,, if J>1. We shall shortly 
compute C, , for both types of boundaries, after which this remaining 
task is easily accomplished by means of relations (41) and (42) for the 
rigid boundary problem and by means of relations (51) and (52) for the 
free boundary problem. 


§ 5. COMPUTING THE EIGENVALUES AS Roots oF A CuBIcC POLYNOMIAL 
The definition (33) implies that if «.8<z<a y,,° then 


ist Wigers i Z ry Zz 
n=1 te (ON =o") & y"(z—a 7°) & v417(2—%& 41°) 
—2 YS <if(z)< YS Ht 
n=N+2 One '( On — Sey) SPOS nad ge (ee ens & y*(z—a 7°) 
cs pe 1 


PERT PERE Ie ER (53) 
ty 41(2—-%y 41°) n=N+2 %n(%pn°—% y°) 
Since all three functions compared in (53) are monotonic decreasing with z 
except across their poles, the root of (32) which lies between «6 and 
%y41° 1s bracketed between the middle roots of the two cubic polynomial 
equations generated by setting the two extreme terms in (53) equal to 
[q(2/-+-3)]-1. The distance between these two middle roots of the two 
cubics increases like N* as N increases. Since 7, y increases like N°, the 
estimate of r, y derived from (53) has a per cent error which decreases 
like N~? as N increases. 

However, computation of the coefficients of the two cubics generated 
from (53) is quite laborious, and 7, y can be computed more simply if / 
is not too large and n=1 or 2. It is necessary to observe that 


~~ sa) ee 
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where 2 
0<0< (55) 
Le a bo 
The two lowest roots 7, and 7, of the cubic 
Zz Ts zZ ( x 1 IK 1 (56) 
a y*(2—ay®) © a9%(Z—0198) nas An®/ 9(20-++8) : 
are then estimates for 7, and 7,, with errors of the order of 
r p 1 eh \e : 
PAisaisy © 


This method of computation was used for C,, and C_,, in tables | and 2 


Ci,3 was computed by Newton’s method. 


Chandrasekhar’s values tor 


C,,,, given in tables 1 and 2 of I, and obtained by his variational method, 
are written to the left of the corresponding exact values in tables 1 and 2 


of this paper. 


Table 1. 


Free Boundary C;,,, 


Variational 1 Exact 2 3 
/ value value 
il 3091-4 3091-180-.0-014 98079 +3 809150 4-50 
‘s 2 5224-4 i 5224-112 40-002 Bi 89189-0 +2-5 
wy 3 8775-0 vi 8774-517 0-004 102124 “44 
+ 13982 13981-34 0-01 124455 +3 
5 21204 21203-58 +0-02 tc 154315 LS 
6 30848 r 30846 891 4. 0-026 
Table 2. Rigid Boundary C,,,, 
n 
Variational 1 Exact 2 3 
ae value value 
1 8047-1 7 8040-07 +0-04 158000 + 350 “1120150 +450 
res 10403 10388-9 +0-2 132500 +350 
3 15132 15104-9 +.0-3 2 143500 500° b. 
io 4 22006 21959 tice 168000 + 600 
5 2) 31315 31238 +2 “5 200790 + 800 oe, 


6 43459 


43342-4+1-9 
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From these tables it is clear that the eigenvalues C,, are too high to 
appear before all C, , convection patterns have begun for / at least up to 6. 
From Chandrasekhar’s tables 1 and 2 it appears that, if the boundary is 
free, convection patterns C,, for / up to 11 set in before any C, » begins ; 
and, if the boundary is rigid, the C, , patterns for / up to 8 set in before any 
C,, pattern appears. It should be noted, however, that an attempt to 
use tables 1 and 2 to predict when various patterns of convection begin 
meets with the difficulty that the linearized Navier-Stokes equation no 
longer describes the flow once any convection pattern has acquired a 
finite amplitude. 

From the tables 1 and 2 observe that Chandrasekhar’s variational 
values for C,,, for / from 1 through 6 are accurate to one part in 10* in 
the free boundary problem and to three parts in 10° in the rigid boundary 
problem. Note also that oh reer 


§ 6. COMPUTING THE EIGENVALUES WITH ARBITRARY ACCURACY 


Very accurate values for the roots of (32) can be computed rather 
quickly by Euler’s method for computing the zeros of J,(z) (Watson 1944, 
p- 500). The advantage of this technique is that, like Newton’s method, 
it can be carried to unlimited accuracy, while less computation is required 
than in Newton’s method. 


a 
Since the infinite product JJ (1—z/x, ,,°) converges absolutely and, in 
n=1 


any bounded region of the complex plane, uniformly to an analytic function 
whose only zeros are simple zeros at z=«.,,°, we can define 


nm? 


w(e)=| 1 garry | (1-35). . (68) 


The function g(z) is analytic in the whole complex plane, g(0)=1, and 
the roots 7,,, of (32) are the zeros of g(z). Using the three facts that 
| 1 41—-%, » | > 77, &, > 1-+4, and that for each / there is a real number ), 
such that for sufficiently large n m(n+31+A,)<a,,<a(n+4l+A,4+4) 
(Watson 1944, pp. 506, 485 and_492 respectively) it can be shown that if 


V gi (2) | i 
Mi as0n i Senn Ree ees) 
l, i { be | 2 |=, (59) 
then lim My, =0. eae ee 


A slight modification of the theorem on p. 137 of Whittaker and Watson 
(1952) then shows that 


n(z)= IT (i=), sata; ‘sn Cee 
n= Tien 
Therefore, if <0 
gizj= 2 gy 42", fo Sails ek oe 


ee 
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it follows that 


. 1 
Ii, z= (—1)* a 


———————-,. ... (63) 
M<N2<... <Np ip nyt, Ly iat set 1). Ne 

and Graeffe’s technique of root squaring may be applied to the power 
series for g(z). The procedure for computing zeros of (58) is to compute 
the coefficients in the power series for g(z) from (58) and (33), and then 
to compute (Graeffe 1837, p. 17 ff.) the power series for the functions 
whose roots are respectively the squares, fourth powers, eighth powers, 
etc., of the roots of g(z). In this connection it is necessary and not difficult 
to show that the power series obtained by applying Graeffe’s algorithm 
to g(z) is simply the infinite product in (61) with r,,, replaced by r,,,?. If 


ice) 
q°@)= ft Oe ee eee A es (64) 


is the power series obtained from (62) by squaring its roots n times, then 


oc 
: _9t a 
Oe ee (ip) ee tae ea es (OD) 


m= 1 
aw 8! 
Since Oar yi ig... aL py = 
, ~m=1 ‘l,m 
converges, therefore if n is large enough it is approximately true that 
nm 
Gr eh cee ge Pose eg =, (66) 
For /=1, the lowest root 7, , of (58) can be computed to eight places by 
taking n=2; that is, by two root-squaring processes. The results are 
C, += 3091-1927--0-0001 for a free boundary 


and C, += 8040-0469+-0-0003 for a rigid boundary. 


§ 7. ConvECTION IN 4 HorizonTaL LAYER oF FLUID 


The method described in this paper can be applied to the classical 
Rayleigh-Benard convection problem in a horizontal layer of viscous 


- fluid heated below (Lord Rayleigh 1916). Supposing that the two boundaries 


of the layer are at z=+4, the eigenvalue problem to be solved consists 
of the equation 
ds 2)" W AW (67) 
te =a ) =—_— oe | op 0 ae no tS 
together with the boundary conditions 
W=0 at z=+4, Poi ee (08) 
2 2 
(Fs -«) W=-0 at z= 4, San te 09) 
and | W‘=0 at a free boundary, | (70) 


or W''=0 at a rigid boundary. | 
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Here A/a? is the Rayleigh number, and ‘a’ is the dimensionless horizontal 
wave number in eqn. (34) of Chandrasekhar (1953 a). For each positive 
number a, there is a critical Rayleigh number, A(a), at which convection 
will set in whose horizontal wave number will be a. The minimum value 
of (a) as a varies from 0 to 00 is the critical Rayleigh number at which 
convection will begin in the horizontal layer. For each a, \(a) is determined 
as the eigenvalue in (67). To compute A(a), let 


RB d* ae 
Define 
[y= (2n+ 1), V, = 20Tr » eat Se meee 
and expand F in Fourier series on [—3, }]: 
F= 2 A, cosp,z+ 2 B, sny2 2) ee 
n=0 n=1 


Then integrating (71) gives 


2 A 
We — (uta-atye seat = fz ae op sin v,z+A coshaz 
+ Bzsinh az+C sinh az+ Dz cosh az. . eras 
But (67), expressed in terms of F, is 
d? 
(zs —at)F= —AW. ~ wae one ee 


Substituting the expression for W given by (74) in (75) and integrating, 
one obtains 


-5=-2 A oe ees ate? 
d Feira Merry ode ce Ga G- =m) 


B C D D 
xzsinh az + an 2" cosh az (5 — im)? cosh az+- i 2* sinh az 
+E sinh az+-G cosh az. oA 3 Ae) oe ee 


Boundary conditions (68) and (69) imply that there exist constants P 
and Q such that 


A=P sinh da, 
B=—2P cosh ja, 


Q— (‘ —a sinh? $a—2 sinh $a cosh =) p 
8a? cosh da ; 
and 
= cosh $a, 


poe (a--a cosh® $a-+-2 sinh $a cosh 4a) 
8a? sinh $a @. 


D=— 2Q sinh ta, | 
| > ae 
J 
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Inserting (73) on the left of (76) and taking Fourier coefficients, 


Ce 2 ake a (= 1) an 
vel ed | 16 cosh? (a) | siete A (79) 
ie : a Kee yee aA 
and B= | 6 sinh? (om) | ; c ci (80) 


Now boundary conditions (70) can be applied at the two boundaries to 
determine \ and hence A, and B,. For example, if both boundaries are 
rigid, differentiating (74) and setting the result equal to 0 at z= +} gives 


i (co) wey n 
ys by»A 4) 1) 


a a+sinha n=0 (Mn? a")? ‘ Sh) 
i! an ae 5} (—1)” 
and —— nn 
se Q o—sinha ,.7 (¥,?+a7)* (82) 
If (79) and (80) be substituted in (81) and (82) the result is 
a+sinha a9 Nu? } 
- Ee cosh? La ‘le in ae aa |=? . (83) 
and 
a—sinh a 2 dv,2 
0 fe sinh? ta el aan |=? . (84) 


If P40, (83) is an equation for A whose solutions, by means of (79), 
produce even eigenfunctions, so that Q=0. If Q0, (84) gives A, which 
produces odd eigenfunctions, so that P=0. 

The case where one boundary is rigid and the other free does not give 
rise to eigenfunctions with a definite parity, so the characteristic equation 
for X, analogous to (83) and (84), is more complicated, although the 
principle is the same. Finally, if both boundaries are free, Rayleigh’s 
original solution (Rayleigh 1916) results from the above analysis. 

Tt would not be useful to reproduce here the eigenvalues one actually 
obtains from the above computations, since the classical Rayleigh problem 
for both rigid and free boundaries has already been solved by a variety 
of methods, and the results are available in the literature. See, forexample, 
Low (1929) or Pellew and Southwell (1940). 


§ 8. MATHEMATICAL SUMMARY 


Because there appears to be a class of problems amenable to the attack 
used in this paper, we shall describe the method in general terms. The 
difficulties produced by singular points will be ignored in this discussion. 

Let B be a linear differential operator of order d with continuous 
coefficients on the interval 0<#<l. Let @B be a class of d linear 
homogeneous boundary conditions on a function and its first d—1 
derivatives at z=0 and «=1. The eigenvalue problem to be solved 
for eigenvalues A and eigenfunctions W is given by the two requirements 
that 

Nee ee (8D) 
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and that W satisfies the boundary conditions B. . oo Sar 


The solution consists of three steps. The first is to find a set of functions 
¢, such that any continuous function W satisfying the conditions ® can 
be expanded in a series 


W= SAd,  . > 2. 


n= 


with constant coefficients 4, which is uniformly convergent on [0, 1]. 
The second step consists in substituting (87) on the right of (85) and 
integrating (85) term by term. The result of the substitution is 


DBW=r J A,dp- 2. ot rs 
n=1 


To perform the integration one must first find the functions (x) which 
satisfy the two conditions 
Dis, (x)=¢,(x) on O<ae<l :, See 
and 
ys, and its first d—1 derivatives are zero at ~=0.* . (90) 


Then one must use the fact that because (87) is uniformly convergent, 
(88) has the particular solution 


W= BA NA, ab, (x . ie aoe 


n 


where the series is uniformly convergent on [0,1]. This can be shown 
by a slight modification of the argument on p. 72 of Ince (1944). Finally, 
one must compute the general solution Q of 


BY=0. s. ns' 10 intl pete at Aa 


This solution Q will have the form 


ad 
Oa)= a Q 34 ;(x) Sse et Sa ae a ae 


where the q; are d linearly independent solutions of (92) and the Q; are 
arbitrary constants. Then the result of step 2 is the equation 


E Agta—Ma)=Q. ss ss (98) 


n=1 


Step 3 requires finding the dual basis for the functions ¢,,—Dv,,, that is, 
the set of functions f(a, A) such that 


[ (ale) Abels 2) dey sss (95) 


* This condition can, in many particular cases, be replaced by less restrictive 
conditions which still ensure the uniform convergence of (91). Clearly 0 ean 
be replaced here by any point in [0, 1]. ’ 
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If k,,,(A) is defined as 


“1 
big A=} Cult Aq; (x) da, Me, eee OG) 


then multiplying (94) on both sides by ¢,, and integrating from 0 to | 
gives 
d 
A= 2 Qjljm(A). Wn ee ee. OL) 
j=1 
If these expressions (97) for A, are substituted in (87), and then the d 
boundary conditions ® are applied to W, one obtains d linear homogeneous 
equations in the constants Q;, the coefficients of which are infinite series 
involving A. Setting the d by d determinant of this system equal to 0 
produces a characteristic equation for the eigenvalues 4. Any solution 
of this equation can be substituted in (97) to give the A,, and hence, 
via (87), the eigenfunction W appropriate to the given eigenvalue. 
The methods by which these three steps can be carried out most con- 
veniently will depend on the nature of the problem at hand. In both 
problems considered in this paper, B is of the form 


D—D,DB,...B a OT re (98) 


.B, 
where 3, is a linear differential operator of order d; ; and the set B of 
boundary conditions contains a subclass € of d, boundary conditions on 
a function and its first d,—1 derivatives. The class 4, can be taken as 
the class of eigenfunctions of the eigenvalue problem 


B,d=2¢ . : . . . . . (99) 
and. satisfies boundary conditions €. fee (100) 


Then the uniform convergence of (87), at least for the second order operators 
3, considered in this paper, is assured by the fact that W must satisfy 
boundary conditions €. Furthermore, in both the problems of this paper, 
there are numbers «, ; Such that ; 


Bib, =%n,iPn> Re Ds Le ou onde ts (101) 


so that we can take 


{= ——- pe ee ee eh LOZ) 


“10.) none An, Dp 


Finally, in both problems considered in this paper, (99) and (100) form 
a self-adjoint boundary value problem so that we can take ¢,, to be ¢,, 
except for a normalization constant. is 
In case the problem under consideration has all the characteristics 
described in the preceding paragraph except the self-adjointness of (99) 
and (100), one can still take y;,, in the form (102), but ¢,, is obtained as an 
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eigenfunction of the eigenvalue problem adjoint to (99) and (100) (Ince 
1944, pp. 210-212). 

In general, if steps one and two can be carried out, step three can always 
be performed in principle as follows* : let @,,() be an orthonormal basis 
constructed from the ¢,—As,, by the Gramm—Schmidt process, so that 


one has 


A(x, A= 2) Cij(A)(h;—Adbs;) . . . . (103) 
j=1 
where ¢,,(A)=0 if i<j. Take €,= 2 ¢,,,(A)O,. Then ¢, is the required 
m=1 


dual basis. 


§ 9. CONCLUSION 


The principal result of this paper can be briefly described. Suppose 
one seeks the eigenfunctions W and eigenvalues A of an eigenvalue problem 
of high order d, of the form given by eqns. (85) and (86), where the boundary 
conditions ® are d linear homogeneous equations in W and its first d—1 
derivatives at 0 and 1. Whenever there is a class of functions ¢,, such 
that any W satisfying the boundary conditions ® can be expanded as 
a uniformly convergent series of such functions ¢,, it is possible to obtain 
a characteristic equation for the eigenvalue A in the form of a d by d 
determinant whose entries are infinite series in X. 

This method is of interest only when, as in the problem of the fluid 
sphere heated from within, it gives a characteristic equation which is 
easier to solve than the characteristic equation obtained by straight- 
forward methods. For Rayleigh’s classical problem of convection in a 
horizontal layer of fluid nothing new is obtained, but in the problem of 
the fluid sphere some new information about the higher radial convection 
patterns with a given angular symmetry can be deduced. Also, in the 
latter problem the method presented here can be used as a check on the 
accuracy of Chandrasekhar’s more convenient variational method. It is 
found that the variational method is quite accurate : two approximations 
are sufficient to give three or four significant figures in the eigenvalues, 
depending on the boundary conditions. 
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ABSTRACT 


A calculation is made of the effect of including the retarded potential 
term of quantum electrodynamics when evaluating the van der Waals 
interaction between two atoms. In disagreement with the conclusions of 
Casimir and Polder the interaction energy at large distances is found to 
vary as R-3. 


§ 1. [INTRODUCTION 


Tue retarded interaction of a pair of atoms was the subject of a paper 
by Casimir and Polder (1948) (hereafter C-P). These authors obtained 
the energy of interaction by means of a perturbation method using the 
electronic charge e as a perturbation parameter. The energy is ~e?. 
In the course of their calculation the unavoidable infinities of energy 
make their appearance but are shown to be of no significance since they 
are independent of the separation of the atoms. The standard result, 
obtained by Eisenschitz and London (1930) considering only the electro- 
static interaction, is that the interaction energy is attractive and oc R-®. 
C—P find that this is modified, at large distances, to an attractive energy 
oc k-7, This result has been widely quoted (e.g. Verwey and Overbeek 
1948). 

A closer inspection shows that the C—P result cannot be regarded as 
a direct consequence of the principles of quantum field theory but 
depends on the particular method of approximation used by them. 
Starting from two atoms and the field of radiation in vacuo they first 
obtain the distortion (~e) of the radiation field due to the presence of 
one of the atoms. This distorted field and the second atom in its ground 
state are then used as the zero approximation for a second stage of 
perturbation which, being carried to terms ~e?, yields energies ~e* as 
required. 

There is no necessity for proceeding in this way in which there is an 
artificial asymmetry between the atoms. As an alternative the original 
unperturbed state can be corrected to fourth order in the energy in a 
completely symmetric way by systematic application of the Rayleigh— 
Schrédinger perturbation theory. 


* Communicated by the Author, 
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‘It is, moreover, no foregone conclusion that the two methods should 
yield the same result. It is well known that the solution of differential 
equations by the method of successive approximations, in particular by 
expansion in powers of a small parameter, does not lead in general to 
unique results. An instructive example was recently given by Kronig 
and Brujsten (1950) who showed that an apparently innocuous problem 
in thermal conduction may become insoluble if the zero approximation 
is not carefully chosen. 

In the present problem it is quite likely that the two methods should 
not agree. The zero approximation of the conventional perturbation 
theory is a non-degenerate state. According to the general formulae of 
the Schrédinger theory (see for example Condon and Shortley 1951, 
pp. 30-35) there can be no resonance denominators in the matrix elements 
in any degree of approximation. In contrast the matrix elements given 
by C—P show resonances due to the degeneracy of the distorted field of 
radiation. 

In this paper the conventional perturbation theory for the above system 
has been carried through with a result differing markedly from that of C—P. 
Since these authors assert that they have checked their unconventional 
method against the ordinary theory, the intermediate steps in the present 
calculation are presented in some detail so that there can be no doubt 
that this particular assertion of C—P is wrong. 

Instead of the reduction of the attraction at large distances as found 
by C—P it is now found to be increased by the retardation, although the 
increase is likely to be insufficient for direct observation. The indirect 
consequences of the increase are, however, by no means insignificant and 
are the subject of a further investigation now in progress. 

The question remains whether the present conventional perturbation 
method or the asymmetric approach of C—P is the correct one. Subject 
to the reservation that, at the present stage of quantum field theory, a 
question of this kind may be premature, it cannot reasonably be doubted 
that the present approach is to be preferred because of the lack of 
singularities in the intermediate stages. 


§ 2. MerHop OF CALCULATION 
The unperturbed state of the system consists of atoms A and B and the 
radiation field all independently occupying their ground states and is 
represented by the wave function 


BRE TR AUIRCCDS 3 S55. A A heal eee Ow) 
The perturbation is represented by the addition to the Hamiltonian of 
the term 2 se ; 
as yy ete: . ; 
H'=Q+ 3 ( <p At soa!) (2) 
eds q, -Rq,-R 
where O Lae = See eee eo. (3) 


4U2 
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is the electrostatic interaction term whose effect is determined in the 
standard van der Waals calculation 


imck,q,=—eap. sa pan toe eer 
A 
denotes summation over all the electrons of atom A. 
A 
hek,=E,—Ey. o 4 Lee Se ee 
E,, E, are the ground and excited state energies of an atom. 
A= 2'c0,e(k,A){Ag, exp (tk. r)+ Ay, exp (—tk.r)}.  . (6) 
kA 


Ay, and Aj, are respectively the annihilation and creation operators of 
quantum electrodynamics with the properties 
Ag al(. Minas + SHV (NM... (N—-Dea- +») | 
Ag... Mga JH=V(N+IE(. (N+ a+ + +) 
C,=(27h/wL)'? (it having been assumed for normalization purposes 
that the system is enclosed in a cubical box of side L where L— o finally). 
e(k,A) is the unit vector representing the polarization of the photon 
whose wave vector is k. A denotes the two polarization possibilities 
corresponding to each value of k. 
It is convenient to regard the electronic charge e as the perturbation 
parameter and the full perturbed Hamiltonian may be written 


(7) 


A =8,+Hf' =8,+-eH,-@i . 50 eee 
where eH ,=— Perie oy Ve ee 
CT a eee mat A alse en 


(N.B. Q contains the factor e?.) 
The perturbed ground state wave function is assumed to be of the form 


P = Pot e 2 Ast ap sPyp ae PLD ih ewe + Ese S S (1 1 ) 
n=(r, 8, t) n 
and the perturbed ground state energy may similarly be expressed in the 
form of a series Hy’ =H, teh, ete, 0 Se ee 
where it may be shown that the first four perturbed terms are given by 
Et) Eh nae . . . . . ° . . . . . . . . . (13) 
ne aera amp (alent (14) 
Ky 3) — Xl » Ho B, (A don not (H von(H 2)nos 


5 2 FLEE EF aloa Haan Hadno 


n#£0 n' +0 ( 
1 
7 (1 s)oo ~ iE, (Avon. 
Hy) =(a)-+()+-(c)+-(d)-++(2), © 9 2 nt ee 


(15) 
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where 


1 
(Raps 
n#0 n’' 40 (E,—E,, E,—#,,) 


) 
x {(H5)on(H1)n n(Hy)not (H4)on(H 2), nH) no (Hy)on(H1)nn(H2)no}> 
p>, 


n#z9 


I 
= — {ada 2 aaa, Who PY 2 camel Heal 


; i 
(¢) E,—E, (eheon (2, 
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(These formulae include the effects of renormalization of the wave 
function and agree with those quoted in Condon and Shortley 1948.) 

The following simplifying assumptions (identical with those of C-P) are 
made : 

(a) Atoms A and B are identical and so far apart that exchange effects 
can be neglected. 

(b) The ground states of the atoms are S states (in other cases there is 
little interest in the van der Waals forces). 

(c) The only effective excited state is the first accessible triply degenerate 
P state (all other states will contribute only to a lower order of magnitude). 

(d) The factor exp (tk . r) can be considered constant for the integration 
over the coordinate space of either atom (this is a common procedure in 
radiation theory and neglects quadrupole and higher order moments). 

With these assumptions no difficulties are encountered in the evaluation 
of the various matrix elements which are found to have the following 


values : 
(1 1)o9=9; ee eee 
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— —icd,{k, exp (—tk- Ry)(4)r08s0+ Xs XP (—tk. R3)(4) 09,05 : e(k,A) 
(n=(r; 8; k,A) n’=(0; 9; 0)) 
=ic0 Aky exp (—ik.Ry)(q)ordoy+ hs’ exp (—ik g Ry)(Q)ovFor} ‘ e(k.A) 
(n=(0; 0; k,A) n'=(r’ ; s’ ; 0)) ete. 
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e2 
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(N.B. n=(r; 8; k,A-+k’,A’) in (22) or (23) and a further factor 1/,/2 
is required if k,A =k’,X’.) 
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The perturbation energy (#)’—,) contains a number of terms arising 
from the separate interactions between atom A and the field and between 
atom B and the field. These are omitted as being irrelevant to the problem 
of the interaction between the atoms. The interaction energy calculated 
by modifying the perturbation energy in this way will be termed JZ. 

Following C—P the coordinate axes are so chosen that the axis of 
separation of A and B lies in the z direction and the degenerate p wave 
functions are so selected that the corresponding dipole moments lie in 
the «, y and z directions respectively. In the course of the calculations 
summation over k is replaced by an integration over k space with the 
correspondence 


ogee ne 
d et dk 
The equalities Sse hy es (oe F= i l(q)o,[? (26) 
2mc® ,. 5 on ae Pt aN? oe 
and dQ€ P(k;Aj= 1 —k?/k? . . . . . . . (27) 


A 


are utilized and it is convenient to make the substitutions v—kR and 
Gah. 
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§ 3. RESULTS 


On examining the perturbation terms in detail it is readily seen that 
E,'Y, Ey and (e) of EZ, are all zero. In addition 1) is eliminated 
by the removal of the separate (A—F) and (B—F) interaction terms. 
Of the terms to fourth order this leaves only (a), (b), (¢) and (d) of Hy. 
These may be divided into three groups: those arising from the purely 
electrostatic interaction, those from the purely retarded field interaction 
and the mixed terms. 

There is only one representative in the first group ; it comes from one 
half of term (c) and is given by 


(28) 


This is, of course, identical with the Eisenschitz—London result. 
The mixed terms arise as part of (b) and combine to give 


4q?2 (” il a : 3cosx 3sinx 
Ky Xiz | de ap) (sine ee “) » (29) 


and this expression is C-P’s 43; (given in their eqn. (40) in different 
form). 

Since X, and X, are identical with the corresponding terms of C-P 
there will be no further discussion concerning them. It is in terms arising 
from the pure retarded field interaction that differences arise. These 
will be considered in greater detail. 

Term (a) can derive from a variety of schemes of virtual states ; one 
such scheme is 


el 
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and the corresponding contribution to AE is 
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From the method of choosing the atomic wave functions this 
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and replacing the summations by integrations over the spaces of k and k’ 
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when 7 represents the x or y direction and 
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RR OAR? 
when it refers to the z direction. 
Written in terms of the dimensionless quantities x, x’ and a and allowing 
for the identicality of the atoms the expression becomes 
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An identical contribution comes from a similar scheme with A and B 
interchanged. The two may be combined and expressed in the form 
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term (b) gives —X 
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Term (d) contains two contributions which cancel exactly with (38) and 
(41). Combining the remaining expressions the net contribution due to 
the pure retarded field interaction amounts to 


and (c) gives 
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for large values of a 


8a? ( dx da = 
X3~ X15 M’'+2NN' ~0-22a9X,.-. . (46 
ve Xisa| | (MM'+2NN") (oo 0-220°X, (44) 
Since Xo~—2X, and 4H=X,4+X,+X%, . . =. - (46) 
this implies that (again for large a) 
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1 
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(In this connection ‘large R’ means greater than a few wavelengths of 
the radiation emitted by the atom in making a transition from the excited 
to the ground state. This will normally be of the order of 10~° cm.) 

It will be noted that there are no infinities in the final result. Their 
disappearance is a consequence of ignoring contributions from the separate 
(A—F) and (B—F) interactions. C—P neglected the infinities on the 
grounds that they are independent of R; this is legitimate reasoning 
but the present procedure brings out the more fundamental aspect that 
there are none of them specific to the interaction between the atoms. 

The present result is to be contrasted with the conclusion of C-P who 
found that 4H varies as (—1/R*) at large R. It should be emphasized, 
however, that in each case the short range approximation agrees exactly 
with the Kisenschitz—London result. 

An obvious extension of the problem is to determine the interaction 
energy of two atoms moving with constant momentum. This is 
under consideration and it is hoped to publish the results in a further 
communication. 
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ABSTRACT 


A search for the existence of K-capture in 3°Cl has been made using 
a proportional tube spectrometer with labelled methyl chloride as a 
gaseous source. A peak which was observed in the beta spectrum 
corresponded to the K-radiation of sulphur and indicated the presence 
of K-capture. The branching ratio K/B- was found to be (1:7-+0-1)%. 
A scintillation spectrometer failed to detect any gamma-ray lines above 
20 kev with intensity greater than 10-* of the beta intensity, so the 
observed K-capture transition goes to the ground state of °°S. 


§ 1. INTRODUCTION 


THE isotope °®Cl is known to be an electron emitter with a half-life of 
(3-08-L0-03) x 10° years (Bartholomew, Boyd, Brown, Hawkings, Louns- 
bury and Merritt 1955) and with a maximum energy of (714-15) kev 
(Feldman and Wu 1952). No gamma-rays with intensity greater than 
5% of the beta particles are present (Wu, Townes and Feldman 1949). 
The spins of the initial and final nuclei are 2 and 0 respectively (Townes 
and Aamodt 1949) and the ft-value is 3x10 sec. The transition is 
second forbidden and the spectrum shape has been extensively studied 
by Wu and co-workers (Longmire, Wu and Townes 1949), as it allows 
a detailed analysis of the kind of interaction (see Wu 1955). The 
observed data can be fitted on the basis of a no parity change transition 
and a tensor interaction. 

36C] is the lightest nucleus which has two neighbouring stable isobars, 
and it might decay both to 3%A by electron emission and to *°S by 
positron emission or electron capture. However, no positrons have 
been observed in experiments which could have detected 10~* positron 
per beta disintegration (Wu, Townes and Feldman 1949). Recent mass 
measurements (Wapstra 1955) give 380 kev for the energy available in 
the transition to 36S. Electron capture is, therefore, possible and a 
search for this was made in the present work. 


§ 2. EXPERIMENTAL PROCEDURE 


Since there are no gamma rays, K-capture can only be detected by 


direct measurement of the sulphur K x-rays, which have an energy 
apa Od Seats iat Dotan hae ORM ENE EPO iy he ee VS te 
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of 2:3 key. For this purpose the proportional tube is the most suitable 
type of spectrometer. However, a solid source would not give reliable 
results since 95%, of the K-capture transitions must result in emission 
of Auger electrons and large errors due to scattering and absorption 
would occur. By using a gaseous source in the proportional counter 
(see Curran 1955) these errors are eliminated, and in the present 
experiments this technique was employed. 

Chlorine itself or hydrogen chloride could not be used as gaseous 
sources since they reacted with the counter walls. However, methyl 
chloride was found to be very suitable. On opening a counter con- 
taining this source to a vessel of equal volume, the activity was exactly 
halved, indicating that none of the source was on the walls. 

The chlorine 36 was supplied by Atomic Energy of Canada Limited, 
as a processed isotope of specific activity 120 u.C/g of chlorine in the form 
of hydrochloric acid free from radioactive impurities. Chlorine liberated 
from the acid was used to prepare methyl chloride by photochlorination 
of methane. Hydrogen chloride and chlorine were removed by reaction 
with sodium bicarbonate and mercury. 

A brass counter of diameter d=10 cm and of sensitive length /=58 em, 
fitted with field correcting tubes, was.employed. It was filled with an 
argon—methane mixture to pressures up to 2 atmospheres. A suitable 
amount of methyl chloride was introduced by diffusion after filling was 
complete. 

The electronic equipment included a Hutchinson—Scarrott maulti- 
channel kicksorter. The 9-2 kev x-rays from an external 74Ge source 
were used for energy calibration. In addition the resolution of the x-ray 
line provided a convenient check of the operation of the counter. No 
change in the shape of the calibration peak was observed when the 
methyl chloride was introduced. 

Since the beta particles from the chlorine produced very large pulses, 
counting rates were kept below 100 counts/sec, and the cosmic-ray 
background was reduced by extensive lead shielding. 

As no change of counting rate was detected in a period of 2 weeks 
the absence of radioactive impurities was established. The possible 
impurities are of relatively short half-life and they would have been 
revealed by a decay of counting rate. 

Part of a typical spectrum, in the range between 0-5 kev and 4 kev 
is shown in fig. 1. A peak at about 2:4 kev, superimposed on a continuous 
spectrum due to the beta particles, is evident. In order to obtain a more 
accurate energy calibration in this region a trace of ?7A was allowed to 
diffuse into the counter after the spectrum measurements had been 
completed. This source decays by K-capture and emits the 2-8 kev 
radiation of chlorine. The 3’A calibration peak obtained is shown dotted 
in fig. 1, the scale on vertical axis being reduced by a factor of 100. This 
shows that the peak from the **Cl source had exactly the same energy 
as sulphur K-radiation. The ratio of the number of counts in the peak 
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to the total number of counts from the source was 1:65°4. In determining 
the area of the °6Cl peak the base of the peak was established by comparison 
with the 3’A calibration spectrum. 


The presence of the sulphur x-ray peak is definite evidence for K-capture 
of 8°Cl since in the above experiment such a peak cannot be produced 


Fig. 1 
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Part of the spectrum of 36Cl showing the peak due to sulphur x-rays. The 
number of pulses in the whole spectrum is 182 000. The dotted curve 
is a calibration peak from an *’A source. 


in any other way. It is worth pointing out that x-rays from excited 
atoms produced in the *°Cl beta decay have an energy equal to that of 
argon K x-rays, and are not detected with a gaseous source since they 
are effectively simultaneous with the beta particles. 


§ 3. BrancHine Ratio 


In order to obtain the branching ratio between K-capture and beta 
emission the value given above for the relative number of counts in the 
peak compared with the number of counts in the total spectrum must be 
corrected for errors introduced by the finite dimensions of the counter, 
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The number of counts in the peak is affected since some of the x-ray 
photons escape from the sensitive volume of the counter, while some of 
those produced outside cause photoelectrons within it and are detected. 
Further, the observed number of beta disintegrations is too large on 
account of electrons entering the sensitive volume of the*counter from 
the ends of the counter. Both these corrections can be calculated 
rigorously, but in the present experiments they are relatively small and 
it is sufficient to estimate them approximately. 

Due to the Auger effect x-ray photons are emitted in only 5-6% of the 
K-capture disintegrations. For the sulphur K x-rays the half-thickness 
in argon at 2 atmospheres is 0-54 cm. Since the distance between the 
boundaries of the sensitive volume and the ends of the counter is h=6-5 em 
which is large compared with the half-thickness, photons entering the 
counting volume compensate for those leaving it. The probability of 
escape of photons through the side walls of the counter has been estimated 
to be 5°%%. The total correction to the K-capture due to the photons is, 
therefore, 0-3°%. 

In the case of the beta particles the error in counting rate is due to 
that part of the source which is outside the sensitive volume of the 
counter, and is evidently less than h/l, which in the counter used is 11%. 
In a first approximation for the calculation of the error the absorption 
of electrons can be neglected since f is appreciably smaller than the 
average range of the electrons, and effective solid angles can be approxi- 
mated for different elements of volume and integrated. The values so 
obtained can be corrected for the absorption of the electrons on the basis 
of the known beta spectrum. Such a calculation gave a value of 4:5% 
for the beta intensity correction. 

After applying both corrections and taking a weighted mean of several 
measurements, the branching ratio K/S~ was found to be (1-7-0-1)%. 

It should be remarked that no measurement of L-capture, which will 
increase the total branching ratio slightly, has yet been made. 


§ 4. Decay SCHEME 

In order to establish the decay scheme a search for gamma-rays that 
might be associated with the K-capture was made using a scintillation 
spectrometer. A *°Cl source of IC activity in the form of hydrogen 
chloride, surrounded by a polythene beta absorber, was mounted in 
front of a Nal(Tl) crystal of length and diameter equal to 5cem. No 
gamma-ray peaks were observed above the continuous bremsstrahlung 
spectrum. Since any gamma peak in the energy region above 20 kev 
with an intensity greater than 1%, of the total bremsstrahlung would 
have been detected, the experiment sets an upper limit for gamma-ray 
intensity of 10-* gamma per beta particles. 

However, any low energy gamma rays present would probably be 
strongly internally converted. The energy region below 50 kev was, 
therefore, investigated with the proportional counter containing the 


K-capture in the Decay of Chlorine 36 1341 


gaseous source. Apart from the K-capture peak no other peaks were 
observed above the beta spectrum. However, since the spectrum is 
completely distorted and is recorded mainly in the region below 100 key 
(as most of the beta particles spend only a small fraction of their energy 
in the sensitive volume) peaks with energy above 15 kev of intensity 
less than 1% of the total counting rate could not have been distinguished 
from the large background. Even if any such converted gamma-rays 
are present the K-capture associated with them would not be recorded 
in the peak shown in fig. | because the counter would give a pulse 
corresponding to the sum of both energies. 

It may be concluded that the observed electron capture transition 
which occurs in 1-7°% of all disintegrations goes to the ground state 
of 36S and the decay scheme is as shown in fig. 2. 


Fig. 2 
36 
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Decay scheme of °6Cl. 


§ 5. REMARKS ON THE USE OF A SOLID SOURCE 


Some results obtained with a very thin solid source of ?°Cl prepared by 
reaction of chlorine gas with the walls of a soft glass counter of external 
cathode type are worth mentioning. As 30g of chlorine were deposited 
on a surface of 1000 cm? the average source thickness is 0-03 g/cm? if 
uniform spreading is assumed. The observed K-capture peak had an 
intensity of less than 10%, of that obtained with a gaseous source in the 
same counter and was shifted to a lower energy. The result shows that 
even very thin solid sources deposited on the wall are not reliable in this 
very low energy region. Further work on this phenomenon will be 
discussed separately. 


ACKNOWLEDGMENTS 


The authors would like to express their thanks to Professor P. I. Dee 
and to Dr. S. C. Curran for helpful discussion and encouragement. One 
of the authors (R. W. P. D.) would like to acknowledge the financial 
support of the Department of Scientific and Industrial Research, while 
the other (A. M.) is indebted to the University of Ljubljana for leave of 
absence, 


1342 On K-capture in the Decay of Chlorine 36 


REFERENCES 


BARTHOLOMEW, RosaLiE M., Boyp, A. W., Brown, F., Hawxinas, R. C., 
Lounssury, M., and Merrirt, W. F., 1955, Canad. J. Phys., 33, 43. 

CuRRAN, 8. C., 1955, Proportional Counter Spectrometry in Beta and Gamma-ray 
Spectroscopy, edited by K. Steqpann (Amsterdam: North-Holland 
Publishing Company). 

FELpMAn, L., and Wu, C. S., 1952, Phys. Rev., 87, 1091. 

Lonemirg, C., Wu, C. S., and Townss, C. H., 1949, Phys. Rev., 76, 695. 

Townes, C. H., and Aamopt, L. C., 1949, Phys. Rev., 76, 692. 

Wapstra, A. H., 1955, Physica, 21, 403. 

Wu, C.8., 1955, Haperiments on the Shapes of Beta Spectra in Beta and Gamma- 
ray Spectroscopy, edited by K. StegBAHN (Amsterdam: North-Holland 
Publishing Company). 

Wt, C.S., Townss, C. H., and FetpMan, L., 1949, Phys. Rev., 76, 693. 


E L343ere| 
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ABSTRACT 


It is shown that a boundary containing dislocations of two systems 
can, in general, form on two, and only two, planes, which have simple 
indices, and that the rotation axes are simple crystallographic directions. 
Most of the experimental evidence on misorientations in deformed crystals 
can be explained in terms of such simple boundaries. 


§ 1. INTRODUCTION 


ReEcENT work (Hirsch and Kellar 1952, Gay and Kelly 1953 a, b) has 
shown that many of the dislocations left inside a metal after deformation 
are arranged along surfaces forming low angle boundaries between 
neighbouring regions of crystal, these latter being called subgrains, cells 
or particles. Such low angle boundaries are observed after either high 
temperature or room temperature deformation. The effect of these 
particles on the mechanical properties is not clear, although various 
suggestions have been made (Bragg 1949, Mott 1952, Gay, Hirsch and 
Kelly 1954). To understand their effect, the nature and method of 
formation of the boundaries must be known. While the size of the 
particles is generally well established, little is known about the mis- 
orientations and the arrangement of the dislocations in the boundaries. 
In the course of the deformation of single or polycrystals, slip may take 
place on one or more slip systems in the same region of crystal. It is likely 
that the dislocations will tend to form surface arrays in order to reduce the 
strain energy. In this paper an investigation is made of the types of 
stable boundary that can be formed with dislocations of one or two systems : 
certain special cases involving more than two sets of dislocations will also 
be considered. The experimental data are discussed in the light of the 
theoretical predictions. It is hoped in a later paper to consider the case 
of slip on three or more systems. 


§ 2. THEORY 


Frank (1950) has shown that the dislocation content of a plane low 
angle boundary can be obtained from the equation 


d—(r x 1)2 sin 4a ae Lae er eas Soc SL) 
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where d is the vector sum of the Burgers vectors of the dislocations inter- 
sected by a line r lying in the boundary, | is a unit vector parallel to the 
axis of rotation and « is the angular rotation about this axis. It follows 
from (1) that d is normal to | and r for all r. Further the dislocation line 
must lie in the boundary plane and in the slip plane which is assumed to 
be one of the commonly observed slip planes.* These conditions define 
the boundary plane and the axis of rotation of a boundary formed with 
dislocations of one or two systems. The actual dislocation content may 
then be found easily by application of eqn. (1).T 


Case 1. Slip on One System 


Since r is normal to d for all r, and d is fixed in direction, the boundary 
plane must be normal to this direction, i.e. the direction of the Burgers 
vector, and so is normal to the slip plane. The dislocation lines are all 
parallel to the intersection of the slip and boundary planes: if we take r 
parallel to the dislocation lines, d=0, and so | must also be parallel to this 
direction. 

This boundary is a pure tilt boundary. 


Case 2. Slip on Two or More Systems in One Plane 


As in the previous case, if the boundary is not parallel to the slip plane, 
it must be normal to it. [is then again parallel to the dislocation lines, 
and the orientation of the boundary is determined by the relative numbers 
of dislocations of the different types, which can vary. The boundary is 
again a pure tilt boundary. 

If the boundary is parallel to the slip plane, since r and d are both 
always in the slip plane, | is normal to this plane. The boundary is a pure 
twist boundary. 


Case 3. Slip on Two or More Systems with One Burgers Vector 
ris normal to d and, as in case 1, d is fixed in direction, so the boundary 
plane must again be normal to the Burgers vector of the dislocations. 
The boundary is a pure tilt boundary, but the axis of rotation may take 
up any position within this plane, depending on the relative contributions 
to the boundary of the dislocations on the various slip planes. 


Case 4. Slip on Two Unrelated Systems 
(a) If the dislocations are not all parallel it is convenient to choose r, 
parallel to one set of dislocations, so that it cuts only the other set of 
dislocations. dj, is then parallel to b,, the Burgers vector of the second 


* For boundaries formed under conditions where an unlimited amount of 
climb is permitted, see Frank (1955). 

} Note that for a given pair of systems, the ratio of the numbers of dis- 
locations of each system is fixed, 
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set. Then it follows from (1) that I is normal to both Burgers vectors. 
The direction of the dislocation lines is defined since r,; must be normal 
to b,: the boundary contains both sets of dislocation lines, and therefore 
it also is uniquely determined. 


(6) If the dislocations are parallel to each other, they must be parallel 
to the intersection of the slip planes and it follows that I is also parallel to 
this direction. Since | lies in the boundary, d must be perpendicular to the 
boundary. Suppose that the boundary contains 7, 2, dislocations per cm 
normal to the dislocation lines with Burgers vectors b,, b, respectively. 
Then (n,b,-+-n,b,) must be normal to I, so nyb,.l=—n,b,.1. Thus 
the relative dislocation content of the boundary is determined, and, 
since (n,b,-+-n.b,) is normal to the boundary, the orientation of the 
boundary is uniquely defined. 


Case 5. Slip on Several Planes with coplanar Burgers Vectors 


~ For all positions of rin the boundary, d is in the plane containing the 
Burgers vectors ; hence | is normal to this plane. No simple gener alization 
seems possible about the orientation of the boundary plane. 

If the slip planes lie in a zone, a special case may arise in which the 
dislocations are all parallel to each other and to the zone axis. The axis 
of rotation is also parallel to this direction. The boundary is a pure tilt 
boundary. 


§ 3. APPLICATIONS 
3.1. Hexagonal Systems 


At room temperature the operative slip systems in hexagonal metals 
are (0001) [1210], while at high temperatures slip can apparently also take 
place on (1010), (1011) [1210] (Boas and Schmid 1930, Kolesnikov 1938, 
Cahn, Bear and Bell 1954). The stable boundaries for various combinations 
of slip systems are given in table 1. The boundary of case 1 corresponds 
to the frequently observed ‘ bend plane ’ (Hess and Barrett 1949, Orowan 
1942). 

nate Bear and Bell (1954) report the presence of pure tilt boundaries 
of arbitrary orientation normal to the basal plane in polycrystals of Zn 
deformed by creep at 250°c, the axis of rotation lying in the basal plane 
(fig. 1). This corresponds to case 2. It should be noted that although 
the boundaries are approximately pure tilt boundaries, in general each 
boundary consists of at least two sets of mixed type dislocation with 
different Burgers vectors. (The description given by the authors is 
incorrect.) It is interesting to observe, however, that the mean orientation 
of the rotation axis appears to be about 5° off the expected (0001) plane. 
For polycrystalline Zn deformed by creep at 350°c, Cahn, Bear and Bell 
report the presence of tilt boundaries with rotation axes shown in fig. 2, 
and the occurrence of slip on (1010) and (1011) planes. The expected 
positions for the rotation axes are also marked on fig. 2. Although a 


4X2 
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number of points lie close to the expected positions, it is clear that some 
of the results do not conform with the theory. 


Table 1. Rotation Axes in Hexagonal Metals 


Slip Slip Rotation Boundary Type of 
plane direction axis plane boundary 
Case 1 (0001) [1210] [1010] (1010) tilt 
Case 2 (0001) several [1210] | (a) [hk.o] [kh.o] tilt 
(6) [0001] [0001] twist 
Case 3 (0001) [1210] [101 7} (1210) tilt 
(1010) [1210] 
(0001) [1210] [101 7] (1210) tilt 
(1011) [1210] 
Case 4 (0001) [1210] [0001] (1120) tilt 
(1100) [1120] 
(0001) several [1210]} [1210] (1010) tilt 
(1010) [1210] 
(0001) several [1210]} [1210] (1010) tilt 
(1011) [1210] 
A (1010) [1210] [01.1] (1010) tilt 
(0111) [2110] 
B (1011) [1210] [11.1] (1100) tilt 
(0111) [2110] 
CG (1011) [1210] [12.1] (2110) tilt 
(1101) [1120] 
Case 5 one or - * 
several several [1210] [OO001] depends on dilslocation 
(0001) content 
(1010) 
(1011) 


Note.—A, B, C are marked on fig. 2. 


The discrepancies between theory and experiment may be explained 
in the following way. It has been assumed that the dislocation lines 
are straight and that they lie completely in their usual slip plane. It is 
possible, however, that the mean direction of a dislocation line may not 
lie in any simple crystallographic plane if it contains a large number of 
ec Figure 3 shows a possible arbitrary configuration of a dislocation 
ine. 


eS eee 
* The jogs might be formed by cross slip or climb, 
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000! 


oro 


Tilt axes in Zn after creep at 250°c (Cahn, Bear and Bell 1954). Shaded circles 
refer to axes deduced from photographs with only two sub-spots to each 


Laue spot. 


8O° FROM 
OOO! 


olio 
Cahn, Bear and Bell 1954). Rotation axes 
Letters refer 
Shaded circles refer to axes deduced 


from photographs with only two sub-spots to each Laue spot. 


Tilt axes in Zn after creep at 350°C ( 
predicted by theory are madicated thus Sse ; 
to corresponding case in table 1. 
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A boundary formed with parallel dislocations of this kind, with the 
same or different Burgers vectors, would have an axis lying at an angle 
of, say, 10° to the ideal orientation. Larger numbers of jogs or possibly 
slip on planes such as (1012) or (1013) could account for the larger 
deviations from the simple theory at the higher temperature. 


Fig. 3 


Possible configuration of dislocation line containing a large number of jogs. 


3.2. Face-Centred Cubic System 


At room temperature the slip systems are {111} [110]. Table 2 gives 


all the possible boundary configurations for slip on one or two systems. 
Apart from the well known kink bands of case 1 (Cahn 1951, Honeycombe 
1950) there is little experimental evidence on misorientations in deformed 


Table 2. Rotation Axes in Face-Centred Cubic Metals 
Slip Slip Rotation | Boundary | Type of 
plane direction axis plane | boundary 
Case 1 (111) [110] [112] (110) tilt 
Case 2 (111) several [110] | (a) in (111) | L(111) tilt 
(6) [111] (111) twist 
Case 3 q (UID) [O11] in (O11) (O11) tilt 
(111) [O11] 
Case 4 (111) [101] (a) [111] (131) mixed 
(111) [O11] (6) [O11] (211) tilt 
f (111) {LOL | (a) [OLO| (O10) twist 
[101] (b) [O11] (100) tilt 
9 (111) [LOL] (a) |L11] (111) mixed 
(111) [110] ») [O11] (O11) tilt 
Case 5 | several (111) | coplanar [110] [111] — mixed in 


Notes.—d, 


general 


primary and cross slip. 


if, primary and conjugate slip in compression. 
primary and conjugate slip in tension. 
a(b), dislocations can react to form Cottrell-Lomer locks. 
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crystals. Barrett and Levenson (1940) report the presence of deformation 
bands parallel to (010) for compressed Al crystals. The primary and 
conjugate systems for this type of deformation would result in just this 
type of boundary. Ball (unpublished work) has recently examined with 
a microbeam a polycrystal of Al extended at 300°c. The Laue photograph 
from a region showing slip on two sets of planes consisted of spotty streaks 
indicating a unique axis of misorientation which was found to. be very 
near to [100]. This again corresponds to the type f boundary of table 2. 
Evans, Layton and Wilman (1951) have reported the existence of the axes 
of rotation shown in table 3 in abraded crystals of copper. The boundary 


Table 3 


Crystal surface | Direction of abrasion Axis of rotation 
(110) [110] } [001] 

(110) [O01] [110] 

(110) [112] [111] 


planes were not determined. Wilman has postulated a new slip mechanism, 
‘rotational slip’, to account for these axes of rotation (Evans, Layton 
and Wilman 1951). It is clear, however, that all these axes of rotation 
are expected to occur if slip takes place on two of the usual translational 
slip systems. In particular, the rotation about [111] could be obtained 
with a number of different pairs of slip systems, contrary to the statement 
made by Evans, Layton and Wilman. It should be noted that although 
the [100] and [111] rotations could correspond to twist boundaries, the 
rotation about [110] occurs at a tilt boundary (see table 2). A determination 
of the boundary planes would therefore serve as a test of the ‘ rotational 
slip ’ hypothesis and of the explanation suggested here. 

Wilman (1951) also reports rotation about [001] in a number of beaten 
and annealed foils of face-centred cubic metals (Cu, Ag, Au, Al, Porro): 
He observes that in many cases the boundary angles are such that there 
is a simple periodic fit between atoms on the two sides of the boundary, 
assuming this to be (001). The theory presented in this paper predicts 
(001) as the boundary plane, so Wilman’s observations may be considered 
as evidence for the existence of cusps in the curve of energy versus boundary 
angle, as predicted by dislocation theory (e.g. Read 1953). 


3.3. Body-Centred Cubic System 

“At room temperature the operative slip systems appear to be (110), 

(112), (123) [111], and of these (110) [111] is probably the most important. 

Table 4 gives the type of stable boundary for slip on (110) [111] systems, 
and a special case for the other systems. 

Deformation bands with planes parallel to (111) 

for Fe and f brass. There is very little evidence on boundar 


have been observed 
ies containing 
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more than one set of dislocations. Barrett and Levenson (1939) report 
the occurrence of deformation bands parallel to (100) and (111) for crystals 
of Fe deformed in tension (cases k, h). Evans, Layton and Wilman (1951) 


Table 4. Rotation Axes in Body-Centred Metals 


Slip Slip Rotation | Boundary | Type of 
plane direction axis plane | boundary 
Case | (110) [111] [112] (111) tilt 
(112) [111] [110] (111) tilt 
(123) [111] [541 } (111) tilt 
Case 2 (110) [111] and [111]] (a) in (110) | L(110) tilt 
(6) [110] (11 twist 
Case 3 h (110) [111] in (111) (111) tilt 
(101) [111] 
Case 4 (110) [111] (a) [011] (111) mixed 
(101) [111] (b) [111] (011) tilt 
- (110) [111] (a) [110] (131) mixed 
J (101) [111] (6) [111] (112) tilt 
k (110) [111] (a) [101] (O01) mixed 
(110) [111] (6) [001] (010) tilt 
Case 5 one or_ mixed in 
several (110) 111] [110] — general 
(112) [111] 
(123) 
Notes.—h, | primary and cross slip. 


7, primary and conjugate slip in compression. 

k, primary and conjugate slip in tension. 

k(b), dislocations may react in a manner analogous to the formation of 
Lomer locks in the face-centred cubic system. 


have observed the rotation axes given in table 5 for abraded crystals of Fe. 
These misorientations may again be readily explained in terms of 
boundaries consisting of two sets of dislocations (table 4). 


Table 5 


Crystal surface | Abrasion direction 


Rotation axis 


(110) [110] [001] 
(110) [112] [111] 
(001) [110] [110] 
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§ 4. Discussion 


It appears that many of the experimental data obtained so far on the 
nature of dislocation boundaries in deformed crystals can be explained 
on the assumption that slip takes place on one or two of the well-established 
slip systems, and that Frank’s formula is applicable. This latter formula, 
which is used as the criterion for stability of the dislocation boundaries 
in this paper, gives the arrangement of dislocations necessary to fit together 
two unstressed regions of crystals misorientated relative to one another. 
It follows that the dislocations are always arranged in such a way that 
their long range stresses cancel out, and such a configuration has a 
minimum potential energy. This interpretation of the stability condition 
may be illustrated by the following example. Suppose in the case of a 
hexagonal close-packed metal slip takes place on several systems in (0001), 
and suppose that an arbitrary boundary is formed. This boundary has 
in general a resultant screw component parallel to the dislocation lines, 
and an edge component in the slip plane perpendicular to the dislocations. 
The long range stresses of an edge component can be annulled by allowing 
the dislocations to lie in a plane normal to (0001); this leaves a total 
screw component, equivalent to a pile up of screws, the long range stresses 
of which could only be cancelled out by slip taking place with a Burgers 
vector not lying in (0001). Since this is not possible the dislocation wall 
is orientated in such a way that the resultant screw component is zero. 
Frank’s formula gives precisely this condition. 

If there are any external stresses acting on the grains, it is possible that 
the dislocation boundaries may be forced out of their equilibrium position. 
A similar effect may occur if a cell structure is formed, for the minimum 
energy dislocation boundaries of such a structure may not be compatible 
with each other. The long range stresses of such boundaries may be 
eliminated by slip on further systems or climb of a suitable amount. 
Deviations from the simple theory can apparently be explained in this 
way. 

With regard to the cellular structure formed during deformation at 
room temperature (Hirsch and Kellar 1952, Gay and Kelly 1953 a, b) 
there is as yet no clear evidence about the nature of the rotations ; from 
the theoretical point of view, a cellular structure may require slip on more 
than two systems and the stability conditions may have to take into 
account the mutual compatibility of the boundaries. It is hoped to 
consider these points in a later paper. 
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ABSTRACT 


Measurements have been made, using polarized light, of the anisotropy 
of the visible absorption of centres induced in quartz by x-ray, gamma 
and neutron irradiation. The ratio of the integrated absorption when 
the electric vector is parallel to the c-axis to that when it is perpendicular, 
is 1-5. An experiment is described which suggests that the optical 
absorption is associated with substitutional aluminium impurity. It is 
shown that a value of 1:3 would be expected from the O’Brien and 
Pryce model of the magnetic resonance centre which is induced under 
similar conditions. 


§ 1. INTRODUCTION 


THE anisotropic absorption of plane polarized light by absorbing centres 
in solids indicates that the wave function of at least one of the states 
involved in the transition has directional properties related to the axes 
of the crystal. That such an effect occurred in the visible bands (A, 
and A,) induced in quartz by irradiation was suggested by the work 
of Ditchburn et al. (1954). They found that the shape of the absorption 
curve and the values of the absorption coefficients depended on the 
orientation of the crystal relative to the incident beam. 

The results of measurements made using plane polarized light, which 
are described in this paper, show that the centre is in fact anisotropic. 
The anisotropy accounts for the directional effect observed by Ditchburn 
et al. and is in reasonable quantitative agreement with that derived 
from the model which O’Brien and Pryce (1954) use to explain the 
paramagnetic resonance results of Griffiths e¢ al. (1954). 


§ 2. MEASUREMENT OF THE POLARIZATION EFFECTS 


The determination of the complete polar diagram, of absorption as 
a function of the direction of the electric vector E, for an induced 
absorbing centre in quartz is complicated by the natural birefringence 
and optical rotation of the crystal. However, for the major effects these 
factors are shown to be negligible or can be taken into account. 

The types of measurement carried out are summarized in table 1. 
Types (i), Gi) and (v) can all be carried out on the same parallel cut 


* Communicated by Professor R. W. Ditchburn. 
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specimen, that is cut so that the c-axis is in the plane of the polished 
surfaces. We have prepared two specimens each having two pairs of 
polished faces in order that all types of measurement could be made 
on the same specimen. This is particularly important when discussing 
the positions of bands since only then is it certain that all the other 
conditions are identical, in particular the specimen and the irradiation. 
X-ray, y-ray and neutron irradiations were used to induce the absorption 
in specimens cut from three crystals of Brazilian quartz. 


Table 1. Types of Measurement 


7 ae Direction of 
| Type Direction of E incident light 
(i) Plane polarized, E || c-axis ris 
(ii) Plane polarized, E 1 c-axis ia 
(iii) Plane polarized, E 1 c-axis ese 
(iv) Unpolarized, but E always L c-axis iene 
(v) Unpolarized, E at all angles to c-axis 1 c-axis 


All the measurements have been made with a Hilger ‘ Uvispek ’ 
spectrophotometer. For the experiments involving plane polarized light 
a polaroid disc, which was attached to the housing of the exit slit of the 
monochromator, could be rotated in the beam. The polarization of the 
light transmitted by the polaroid exceeded 99-8°% except near the extreme 
limits of the visible spectrum. The plane of polarization could be set 
to about +4°. 

In order to eliminate any effect associated with the rotation of the 
plane of polarization by the crystal and the selective response of the 
photocell, the latter has been determined. The polarizer was set so 
that E was vertical. Let J) be the intensity of plane polarized light 
which, for the first measurement, fell on the cell giving a response kJ, 
where k is a constant. An unirradiated quartz specimen, cut so as to 
give a rotation [«(A)] of 80° at 66004, transmitted an intensity 7, when 
placed in the beam. The response of the cell was then k/,R,, giving a 
measured transmission 7';=7',R,, where 7)=I,/I, and R, is the ratio 
of the response of the cell to plane polarized light having E at an 
angle «(A) from the chosen initial direction. This was repeated for 
settings of the polarizer such that the incident E was » .«(A) from the 
vertical, where n is an integer. When n=1 we have 

: Ledisi ee Ms os 
D y= a: = 7 Roya OF Rog an= rs" 


and for 
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From measurements of 7 the values of R,,) were found for the two photo- 
cells for different wavelengths and some of the results are shown in fig. 1. 
The value of 7’) computed from these results, when (n-+1)x(A)=180, 
agreed within 0-5°% with the value calculated from refractive index data. 
The geometry of the photocell surfaces is also indicated in fig. 1. It will 
be seen that the response was greatest when the incident E had a 
component perpendicular to the surface of the cell. 


Fig. 1 
ae P F 


Incident beam | S|ncident beam 


=o) 


i 
oy, 
ey 
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»=5890A 1 =6600A 


The relative response (R, in %, broken line) of the photocells to plane polarized 
light having E at an angle to the vertical P. The geometry of the 
emissive surfaces is shown above the polar diagrams, the L.H. cell 
being normally used from 2100 to 6500 A and the R.H. cell from 6( 100 
to 10000 A. 


NY 


Various type (iv) and (v) measurements have been reported for smoky 
quartz, e.g. Cohen (1954), and the results are similar to the previously 
mentioned type (iv) and (v) measurements of Ditchburn e¢ al. (1954) on 
irradiated quartz. To compare these results quantitatively with the 
type (i), (ii) and (ili) measurements the polarization introduced by the 
monochromator must be determined. In our case the light is partially 
polarized by the six reflections in the monochromator. 

Measurements were made with the polaroid attachment in the beam, 
the response of the photocell being observed as the setting of the polaroid 
was varied. From the known selective response of the cell it was then 
possible to determine the variation of the relative intensity of radiation 


1356 E. W. J. Mitchell and E. G. 8. Paige on the 


resolved into a given plane of polarization, as a function of the angular 
variation of the plane. A typical result is shown in fig. 2. The orientation 
of the axis of the figure relative to the exit slit was found to depend on 
the wavelength, probably owing to the rotation of the plane of polarization 

in two crystalline quartz lenses in the instrument. 


Fig. 2 


P 


100°, 


The polarization figure of the monochromator at 6100 A. The radia scale is 


in % of maximum polarization. 


§ 3. RESULTS 


The full line absorption curves in fig. 3 correspond to measurements 
of types (i) and (ii) in table 1. It will be seen that the absorption is 
greater when E of the incident light is parallel to the c-axis of the crystal 
than when it is perpendicular. In addition the A, band (2-0 ev) is more 
pronounced, being hardly detectable when E is perpendicular to the axis. 
At —195°c the A, peaks (about 2-7 ev) on both curves have moved to 
lower energies and the A, band on the higher absorption curve is less 
well marked. 

In fig. 4 absorption curves are shown for measurements of types (iii), 
(iv) and (v) together with the curves (broken lines) for (i) and (ii), all 
referring to the same specimen. Using the results for the instrumental 
polarization and the selective response of the cell it was found that the 
instrumental effect arising because of the natural rotation was negligible 
for type (iv) measurements. On the other hand for (iii), in which the 
light was deliberately plane polarized, a small correction had to be made ; 
the curve shown in fig. 4 has been corrected. 

The absorption curves corresponding to (ii), (iii) and (iv) are in good 
agreement showing that the absorption of the centre is isotropic for any 
orientation of E perpendicular to the c-axis of the crystal. When the 
. incident light is perpendicular to the c-axis, as in (ii) and (v), some 
ellipticity will be introduced but the effect is small, as may be seen from 
the transmission when the specimen is between crossed polaroids. Further- 
more, the agreement between (ii) and either (iii) or (iv), in which only 
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Fig. 3 
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The anisotropic absorption in quartz at room temperature (full line) and at 
—195°c (broken line). The upper curves are type (i) and the lower 
type (ii) measurements, 


Fig. 4 
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The anisotropic absorption of quartz at room temperature. The broken line 
refers to type (ii) and the dotted line to type (iii) measurements. The 
three crosses are the predicted type (v) results. 
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rotation is present, shows that E remains predominately in the initial 
direction. Using the curves for (i) and (ii) together with the polarization 
figure for the beam from the monochromator it is possible to predict the 
absorption coefficients for a type (v) measurement. Three points deter- 
mined in this way are shown on fig. 4. It will be seen that they lie on 
the curve for the type (v) measurements. 

Apart from those cases in which all types of measurement have been 
made on one specimen the ratio ;.(E||) : (EL) can be found from type (i) 
and (ii) measurements on a parallel cut specimen having one pair of 
polished surfaces. A summary of results is given in table 2. 


Table 2 


; Position (ev) of 
Specimen : BE||/WEL 
and irradiation oa Pas at 2-6 ev 


N.37 x-ray 
N.36 y-ray 
INGE ai 

N87 
N.14 n 
NLS 7 
N.35 n 
N.39 
N.40 7 
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Although measurements have not been made on the C band (5-6 ev) 
with plane polarized light, the results of the type (iv) and (v) measure- 
ments shown in fig. 4 suggest that the absorption in this band is not 
significantly anisotropic. 


§ 4. Discussion 


It is clear from these results that the A centres give rise to anisotropic 
absorption with respect to plane polarized light. Therefore it is necessary 
to know the crystallographic orientation and the instrumental polarizations 
when comparing the results from different specimens and from different 
kinds of irradiation, obtained in different laboratories. 

A simple comparison may be made for type (iv) measurements since, 
in general, this will be independent of instrumental effects. Unfortunately 
the A, band is not observed in this case so that for a detailed comparison 
type (i) measurements must also be made. Since this involves using a 
parallel cut specimen type (ii) measurements can conveniently be made 
at the same time. 

For the specimens which we have measured, the integrated absorption 
is greater when E is parallel to the c-axis than for any perpendicular 
direction, The concentration of centres involved is the same in the two 
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cases so that the transition probability for the centre must be anisotropic, 
being greater in the direction of the c-axis. An indication of the magnitude 
of this effect is given by the values of (E||): (EL) in table 2. Since 
the half-width does not vary significantly these values are equal, to 
within 5%, to the ratio of the sums of the f-values resolved in the two 
directions. 

The Nature of the Centre 


Brown and Thomas (1952) found that synthetic crystals coloured less 
readily than natural crystals and the results of Mitchell and Paige (1954) 
showed that the A, and A, bands were probably associated with impurity 
since the absorption did not increase with prolonged neutron irradiation. 
It was further shown by Ditchburn et al. (1954) that the effect was due 
to aluminium impurity. 

The aluminium which, after ionization, is effective in producing the 
optical centre is probably incorporated substitutionally, as may be seen 
from the following experiments. Firstly, after heating for 12 hours in 
molten aluminium at 700°c a crystal became black throughout indicating 
that aluminium atoms have diffused for distances of the order of milli- 
metres. Secondly, a non-uniformly coloured crystal, in which the visible 
absorption was in well-defined bands approximately parallel to a major 
rhombohedral face, was heated at 950°c for 24 hours and slowly cooled. 
After a second irradiation it was found that these bands were in the 
same positions in the crystal as before the heat treatment. The simplest 
interpretation of these experiments is that in the first we were dealing 
with diffusion via interstitial sites, while in the second the aluminium 
effective in producing colour centres was present substitutionally. 

That the optical A bands and the paramagnetic resonance effects are 
associated with the same centre has been established by the magnetic 
measurements of Dr. Owen and Dr. Bagguley of the Clarendon Laboratory, 
Oxford, on two specimens in which the A, and A, bands have been 
completely removed. Dr. Owen has shown that the aluminium type of 
resonance was absent in a specimen in which the A bands had been 
removed by optical bleaching, and Dr. Bagguley found the same result 
on a specimen which had been radiation-bleached using a large neutron 
dose (Mitchell and Paige 1954). 

It is concluded, therefore, that both optical and magnetic effects are 
associated with substitutional aluminium. This is in agreement with 
the model which has been proposed by O’Brien and Pryce to account 
for the magnetic results. In their model the substituted Al is assumed 
to be Al- so that in the unirradiated state the binding is normal.* On 
irradiation an electron is released from the defect and is trapped either 
at some other impurity (O’Brien and Pryce 1954) or at a defect of the 
kind which can be produced by neutron damage (Ditchburn et al. 1954). 


Berna yk a ee 
*It may be that the small differences in lattice parameter of different 


quartzes observed by Keith (1950) are due to varying concentrations of 
substitutional Al. 
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Wherever the electron is trapped O’Brien and Pryce suggest that the 
hole spends most of its time on an oxygen adjacent to the Al’. ‘The 
possible optical transitions on this model, using the O’Brien and Pryce 
nomenclature, are sp, to p,, or if the non-bonding sp, is not filled between 
the bonding and non-bonding sp, orbitals, or between p, and non-bonding 
Sp-- 

"The origin of the line broadening is not understood, especially as the 
width is not greatly temperature dependent (Ditchburn e¢ al. 1954). 
However, the anisotropy of the integrated absorption due to the centre 
must be associated with the directional property of the p,, p, or p, 
orbitals since the transition is only allowed by virtue of the s-component 
of the other level. In the discussion of the wave functions given by 
O’Brien and Pryce, p, is the component of the p wave function in a 
direction perpendicular to the planes of the Al-—O+-Si triangles, p, is 
the component parallel to the line joining Al- and Si, while p, is the 
component perpendicular to the line joining Al~ and Si. One bond in 
each of these triangles is at 44° 15’ to the c-axis and the other at 66° 15’. 
Using the data given by Wyckoff (1948) it is found that the normals to 
the six possible triangles make angles of 51° 3’ to the c-axis and the 
z directions 58° 40’, corresponding to the 54° between the Al—Si 
directions (y) and the c-axis. 

We have to compare the measured f-value ratio with 27,/2f,. In 
the parallel direction there are six components of dipole moment MW cos @ 
so that 2/,=C . 6M? cos? 6, where C is a constant. In the perpendicular 
direction, if E makes an angle ¢ with the projection in the basal plane 
of the required direction the components of dipole moment are two of 
M sin @ cos ¢, two of M sin @ cos (6+ 60) and two of M sin 6 cos (60—¢). 
Then 

Sf =C . 2M? sin? 6 {cos? d+ cos? (60+ ¢)-+ cos? (60—¢)}, 
which reduces to C. 3M? sin? 6, being independent of ¢. This agrees 
with the experimental result that the absorption was isotropic for all 
directions of E perpendicular to the c-axis. 

The ratio of the integrated f-values, X/,/2f,, is 2 cot? @ giving 1:3 
when @ is 51° 3’, 1-06 for 54° and 0-74 for 58° 40’. The measured value 
is close to 1-5 (see table 2) so that we reject the possibility that the sp, 
is the higher of the non-bonding orbitals. The agreement with the 
model if p,, is uppermost is considered reasonable. In order to obtain 
1-5 from the model 6 would have to be 49° 6’ which is only about 2° 
away from the normals to the triangles as calculated from the crystallo- 
graphy of undistorted quartz. Random variations about the 51° 3’ 
direction would lead to a higher value than 1-3 since d/d@ (cot? 6) 
decreases with increasing 0. 

The anisotropic optical absorption of the A, and A, bands is therefore 
in reasonable agreement with the O’Brien and Pryce model if the p,, level 
is uppermost. However, the origin of the line broadening and the nature 
of the corresponding electron trap remain uncertain. 
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IXLVII. CORRESPONDENCE 


The Mean Lifetime and Some Properties of 0°-Particles : A Correction 
to Previously Published Data 


By D. B. GayTHER 
The Physical Laboratories, The University, Manchester 13 + 


[Received August 19, 1955] 


Tue Manchester cloud chamber group, led by J. A. Newth, have recently 
pointed out a numerical error in table 2 of a paper by the author (Gayther 
1954) on the mean lifetime of the 6°-particle. The value quoted for «* is 
0-41, whereas the correct value is 0-82. This error does not affect the 
- classification of the V°-decays described in the paper but it does affect 
the values of the momenta of the #°-particles and thus their mean lifetime. 
The corrected values of the various quantities depending on e* are given 
in table 1. The explanation of the symbols is given in Gayther (1954). 


Table 1. Corrected Data on 8 4°-particles and 2 Probable 6°-particles 


(1) (2) (3) (4) (5) | (6) (7) (8) (9) 
ze E W t fis Ts 
Event yB (Mev/c) (Mev) ———--—- --_+ 
10-19 sec 
A1208 0-68 330 100 | 2-2 | 0-6 | 4-0 —| — 
B1843 1-20 590 270 | 2:5 | 0-05] 3-0 —|{; — 
B1990 1-10 550 240 | 1:8 | 0-5 | 43 —| — 
D1967 1-30 640 310: | 1-6: |-1-0 7) 4:3 Se Oe 
D3587 ~7-0 ~3400 ~3000 | 4:5 | 0-1 | 0-5 —| — 
D3754 0-73 360 120 | 3-8 | 0-6 | 6-8 | 3:7 | 48 
D5123 1-30 640 310° ET OF ST ee Oe aes 
D5155 2-50 1200 840° | 1579) 02 > (el O10 Misia 
B3602 _ ~780 ~430 | 51 —} —}| —| — 
D319 —: ~710 ~370 1:7 


The value of the mean life, 7, for the 6°-particle, found from the measure- 
ments on the first eight cases is given by 


1/r=(1-8-E0-8) x 101. sec-? 


i.e. +=(0-670) x 10-1 sec. 


+ Now at A.E.R.E., Harwell. 
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This corrected result may be combined with the data of other groups listed 
in table 6 of Gayther (1954) to obtain a weighted mean value (from 34 
ls r= (1-747) x 10-19 sec. 

In a recent paper Gayther and Butler (1955)* made further use of the 
data on 9°-particles given by Gayther (1954). Owing to the error in the 
original value of <* the values of the momenta and kinetic energies of the 
0°-particles given in table 1 of the G.B. paper are incorrect, thus all the 
conclusions on the production processes of 0°-particles need modification. 
The corrected value of the momenta and energies are given above in 
columns (3) and (4) of table 1. Column (5) contains the corrected. 
weighting factors, W ; the derivation and significance of these factors are 
described in § (3) of G.B. 

The differential kinetic-energy spectrum and angular distribution of 
A°-particles were discussed by G.B. The behaviour of these particles is 
similar to that of the ‘ grey ’ protons from stars. The corrected data on 
the differential energy spectrum for 6°-particles is given in table 2, where 
a comparison is made with A-particles. The shape of the angular 
distribution for the §°-particles given in fig. 3 of G.B. is not significantly 
affected by the change in the weighting factors. 


Table 2. The Differential Energy Spectra for #°- and A°-particles 


<50 |50—-100 | 100-400 | 400-700 | 700-1000 | >1000 


A ee : 
Kinetic energ5 mev | Mev | Mev Mev Mev Mev 


Corrected no. of 
§°-particles ce = — 15 
(actual no. ob- 
served 10) 


1 
bo 
tN 


Corrected no. of 

A°-particles .. 28 10 26 8 — — 
(actual no. ob- 
served 28) 


Superficially the 6°-particles appear to be of higher energy than the 
A°-particles. It is necessary, however, to consider the bias effects at 
low energy due to the size of the cloud chamber and the short lifetime of 
the particles. Reference to fig. 1 of the G.B. paper shows that for 
A°-particles in the energy range (50-100) Mev the detection probability 
ranges from 0-45 to 0-50. For 6°-particles in the same energy interval 
the probability range is 0-30 to 0-42. In this particular energy interval 
a corrected number of 10 A° were observed (actual events 5) and no 
case of a 6°. Thus there is some evidence that the differential energy 
spectrum of 9°-particles probably falls below 100 mev, while Gayther and 


* Subsequently referred to as G.B. 
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Butler concluded that for A°-particles it rises. There is some provisional 
evidence that the spectrum of the 6°-particles is similar to that of the 
7*-mesons produced in stars which, according to Camerini e¢ al. (1950), 
probably falls for energies of less than 100 Mev. 

The analysis of Gayther and Butler dealt with a total of 45 V°-particles, 
28 were identified as A°-particles, while 10 were 6°-particles ; the remaining 
7 cases were unidentified. The relative number of /°- and 6°-particles 
at production can be determined using the detection probabilities W. 
The corrected data give 


N jo: N=(90417) : (2749) 


assuming that all the unidentified events were due to A°-particles. If, 
however, it is assumed they were 0°-particles, then 


N jo: No=(72£16) : (44412). 


Taking into account the possibility of detection bias effects discussed in 
G.B., it now seems that more /°- than @°-particles are produced in nuclear 
interactions in lead. 

A method for finding the frequency of production of A®% and 
6°-particles with respect to the number of shower particles was described 
by G.B.; the corrected results are 


7 O/ . rea Fy O/ 
N ej Nat 3 /0 9 N g/g 13 /O* 


The data refer to V-decays in showers containing less than 9 shower 
particles. 

In conclusion it may be pointed out that if more A°- than 6°-particles 
are produced in energetic nuclear interactions, there must be processes 
other than the double production process. 


a +p>A°+69, 


first reported by W. B. Fowler et al. (1953). One possible additional 
source of A°-particles may be the nuclear capture of K~-particles. Barker 
(1953) and De Staebler (1954) have reported evidence for the processes 


K-+p +/A°+7° 
and K-+n>/°+7-. 


If the K--particles are readily captured in this manner, a considerable 
number of A® may be produced unassociated with 6°-particles. 
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CXLVIII. REVIEWS OF BOOKS 


Advanced Mathematics for Engineers. By H. W. Reppick and F. H. Mituer, 
[Pp. xiv-+548.] (New York: John Wiely and Sons, Ine. ; London : 
Chapman and Hall, Ltd.) (3rd edition, 1955.) Price 52s. 


Tue third edition of this excellent text-book can be warmly welcomed. The 
range of topics is extremely wide and physical applications of the basic mathe- 
matics are well illustrated by problems worked in the text. In addition, lists 
of problems are provided for the student to try, with answers given at the 
end of the book. For the third edition, these lists have been extensively 
revised and the total number of problems has been increased. E.N.F. 


An Introduction to Stochastic Processes, with Special Reference to Methods 
and Applications. By M. 8. Bartuerr. [Pp. xiv+312.] (Cambridge : 
University Press.) Price 35s. 


Tuts long-awaited book is most welcome. Statisticians have been familiar for 
some time with Professor Bartlett’s North Carolina lecture notes ; the present 
work has developed out of these notes into an epoch-making general treatise 
which concentrates on describing and applying stochastic processes: the basic 
mathematical theory, while touched on, is left to be developed more fully by 
Mr. J. E. Moyal in a companion volume which should be published before long. 

The older generation of statisticians, in this country at least, paid little atten- 
tion before the war to the proper methods of dealing with data for which in- 
dependence of individual elements could not be assumed ; for example, analysis 
of economic data could hardly be regarded as satisfactory when done by stock 
methods. That generation, because of the age factor, is having difficulty in 
assimilating the new knowledge in its present rapid state of growth, and so will 
welcome a book which summarizes progress to date; the book, of course, is a 
‘must’ for the younger generation of science students. The greater part of it 
takes us through what must be regarded as a necessary standard treatment of 
the various aspects of stochastic processes. Chapter 5 summarizes some of the 
more formal mathematics of limiting stochastic operations, although a full 
treatment is left to Moyal, while towards the end the basic elements of * clas- 
sical’ statistical analysis are applied to the processes, with particular reference 
to the analysis of time series. 

It is impossible in a short notice to give a proper summary of the main part 
of the book, but mention of the problems dealt with in various scientific fields 
will be a sufficient indication of its all-embracing character. In physics the 
author deals with crystal lattices, diffusion, cascade showers, molecular associ- 
ation in statistical mechanics, and the fluctuation theory of molecular concen- 
trations, with an aside on the Monte Carlo method. In engineering the impor- 
tant theory of communication gets a large share of attention. Genetics is 
represented by the study of genetic recombination, the theory of inbreeding, 
natural selection and population genetics generally. There are important 
population problems in regard to fluctuations, and birth and death processes ; 
in biology growth and mutation in bacterial populations are studied, likewise 
contagious distributions ; in the medical field the theory is applied to epl- 
demics. More generally, there are industrial applications in questions such as 
the renewal of articles with effective random lifetimes, sequential analysis, and 
the theory of queues, as applied to the servicing of machines and control of air- 
craft landings. Examples are given of insurance risks, of the response of a 
learner to a stimulus, and of the varying fertility of an area of land. This will 
be enough to show that the book will be of interest to a much wider class of per- 
sons than the mathematical statisticians, and the author is to be 2 te 


on his industry. 
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Chemistry of the Solid State. Edited by W. E. Garner. [Pp. vitit-417.] 
(London: Butterworths Scientific Publications; New York: Academic 
Press Incorporated.) 50s. 


Tus work is a collection of essays, many of them by Professor Garner, his 
staff and former pupils. It would be better described as the chemistry of 
reaction in the solid state, dealing as it does with nucleation, photochemical 
effects and so on. Several of the articles describe the work carried out in 
Professor Garner’s own laboratory. It is a useful and stimulating volume, and 
the reviewer’s chief feeling on looking through it was that some of the articles 
ought to be expanded into a book. Or alternatively that this new subject, 
the chemistry of reactions in the solid state, needs a new book that would 
attempt a synthesis and compare what is known about the compounds discussed 
with similar information about metals. Perhaps one of the authors of this 
volume will attempt it. They are : 


F. C. Frank—Chemistry of Crystal Dislocations : 

F. S. Stone—Lattice Defects in Ionic Crystals ; 

P. W. M. Jacobs and F. C. Tompkins—Action of Light on Solids : 

P. W. M. Jacobs and F. C. Tompkins—Surfaces of Solids ; 

T. J. Gray—Semi-conductivity and Magneto-chemistry of the Solid 
State ; 

W. J. Dunning—Theory of Crystal Nucleation from Vapour, Liquid and 
Solid Systems ; 

P. W. M. Jacobs and F. C. Tompkins—Classification and Theory of Solid 
Reactions ; 

W. E. Garner—Kinetics of Endothermic Solid Reactions ; 

W. E. Garner and L. L. Bireumshaw—Kinetices of Exothermic Solid 
Reactions ; 

C. E. H. Bawn—Decomposition of Organic Solids ; 

A. R. Ubbelohde—Explosion and Detonation in Solids : 

A. J. E. Welch—Solid-Solid Reactions ; 

J. W. Mitchell—Photographie Process ; 

T. B. Grimley—Oxidation of Metals ; and 

F. S. Stone—Electronic Factor in Chemisorption and Catalysis. 


N.F.M. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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